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Triangle  basis  functions  for  non-modal  formulation. 
Pulse  basis  functions  for  non-modal  formulation. 

Basis  functions  for  modal  formulation. 

x 


Chapter  1 


INTRODUCTION 

The  problem  of  electromagnetic  transmission  through  apertures 
in  conducting  planes  of  finite  thickness  has  been  the  object  of  several 
investigations.  This  problem  is  important  because  of  the  need  to  de¬ 
termine  the  effects  of  wall  thickness  on  coupling  through  slots  and 
holes  in  waveguides  and  cavities.  Several  quasi-static  solutions  have 
been  developed,  which  include  expanding  the  field  in  terms  of  the  char¬ 
acteristic  functions  of  Laplace's  equation  [1],  using  asymptotic  ex¬ 
pansions  to  solve  a  set  of  integral  equations  [2],  and  variational 
methods  [3, A].  These  solutions  are  usually  valid  when  the  radius  of 
the  circular  aperture  or  the  width  of  the  slot  is  small  compared  to 
the  wavelength.  The  effects  of  wall  thickness  have  also  been  measured 
for  some  waveguide  coupling  problems  [5],  and  for  a  long  slot  in  a 
thick  screen  [6,7,8]. 

The  slit  (slot  of  infinite  length)  of  rectangular  cross  section 
in  a  thick  conducting  screen  is  the  object  of  most  of  the  theoretical 
investigations.  The  fields  in  the  slit  region  are  usually  expanded  . 
in  terms  of  the  parallel  plate  waveguide  modes,  while  the  fields  in 
the  half  space  regions  have  various  integral  representations.  Enforcing 
continuity  of  the  tangential  field  at  the  aperture  faces,  one  obtains  a 
coupled  set  of  integral  equations.  The  formulation  and  solution  of 
these  equations  make  use  of  a  variety  of  techniques,  such  as  the  use  of 
Weber-Schafheitlin  discontinuous  integrals  [9,10],  Wiener-Hopf  methods 
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and  generalized  matrix  techniques  [11],  Fourier  integral  methods  [12,13], 
Green's  function  methods  [14,15],  and  the  equivalence  principle  and 
moment  methods  [16].  Another  procedure  is  to  expand  the  fields  in  the 
half  space  regions  in  terms  of  the  characteristic  functions  of  the 
Helmholtz  equation  in  elliptic  cylinder  coordinates  [17,18].  A  doubly 
infinite  set  of  linear  algebraic  equations  is  then  obtained  for  deter¬ 
mining  the  coefficients  of  expansion. 

A  slightly  more  general  problem  is  treated  in  [19],  where  the 
thick  ground  screen  is  an  infinite  slab  with  finite  conductivity. 

Analytic  properties  of  finite  Fourier  transforms  are  used  to  reduce 
the  problem  to  the  solution  of  a  single  variable  Fredholm  integral 
equation  of  the  second  kind.  A  slit  of  arbitrary  cross  section  is 
treated  in  [  20]  whereby  the  scattered  electric  field  is  expressed  as 
an  integral  over  the  electric  current  induced  on  the  infinite  conduc¬ 
tors.  By  requiring  the  total  tangential  electric  field  to  vanish  at 
a  finite  number  of  points  on  the  conducting  contours,  a  system  of 
algebraic  equations  is  obtained  which  determines  the  electric  current 
at  these  points. 

Problems  of  this  nature  are  also  of  interest  in  the  area  of 
electromagnetic  compatibility.  If  the  slit  region  is  filled  with  a 
conducting  material,  it  might  represent  a  gasket  between  two  perfectly 
conducting  walls.  A  different  analysis  than  the  above  mentioned  methods 
is  then  possible  using  approximations  and  models  for  the  transmission 
mechanism  [ 21] . 
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This  report  considers  the  problem  of  a  filled  slit  of  arbi¬ 
trary  cross  section  in  a  thick  conducting  screen.  The  formulation 
presented  here  is  basically  a  specialization  of  that  given  in  [22] 
for  a  class  of  aperture  coupling  problems.  The  basis  of  the  method 
rests  on  the  equivalence  principle,  as  it  is  defined  and  used  in 
[23,24],  which  is  used  to  break  the  problem  up  into  three  isolated 
regions.  This  isolation  is  achieved  by  the  use  of  equivalent  elec¬ 
tric  and  magnetic  current  sources,  which  radiate  in  unbounded  space 
and  are  constrained  to  give  the  correct  fields  in  each  region.  An 
operator  equation  is  thus  obtained  for  each  region  in  terms  of  the 
equivalent  sources  as  the  unknowns.  These  operator  equations  are 
solved  simultaneously  via  a  Galerkin  procedure  [25,26]  after  suit¬ 
able  sets  of  expansion  functions  are  chosen  to  represent  the  un¬ 
knowns.  Some  immediate  advantages  of  this  formulation  are  summar¬ 
ized  as  follows: 

1)  The  necessary  equivalent  sources  exist  on  a  finite 
two-dimensional  contour. 

2)  These  equivalent  sources  radiate  into  unbounded, 
homogeneous  space  and  hence  their  fields  have  a 
simple  retarded  potential  representation. 

3)  The  system  of  three  simultaneous  operator  equations 
has  a  unique  solution  at  all  frequencies. 

4)  The  cross  section  of  the  slit  may  be  an  arbitrary 
shape  and  filled  with  lossy  material. 

The  operator  equations  are  derived  in  detail  in  Chapter  2, 
where  an  argument  for  the  uniqueness  of  solution  is  also  given. 


Each  polarization  is  handled  separately  in  Chapters  3  and  4,  The 
term  transverse  electric  as  used  here  means  that  the  electric  field 
lies  in  the  plane  of  the  slit  cross  section,  and  transverse  magnetic 
means  that  the  magnetic  field  lies  in  the  plane  of  the  slit  cross 
section. 

Once  the  tangential  fields  over  the  aperture  faces  are 
known,  several  transmission  characteristics  may  be  computed.  These 
are  defined  in  Chapter  5,  where  they  are  written  in  terms  of  a 
measurement  vector  and  as  such  are  referred  to  as  measured  quantities. 

A  slit  impedance  and  polarizability  are  also  defined  which,  when  the 
slit  width  is  electrically  small,  characterize  the  slit. 

An  alternate  method  of  solution  is  presented  in  Chapter  6 
when  the  slit  cross  section  is  composed  of  a  number  of  rectangular 
regions.  These  are  viewed  as  a  sequence  of  two-dimensional,  infinitely 
long,  rectangular  cavities,  each  coupled  to  the  other  by  an  aperture. 
The  method  is  an  extension  to  that  presented  in  [16],  and  is  referred 
to  here  as  a  modal  solution.  It  has  the  advantage  that  each  rectangu¬ 
lar  region  may  be  filled  with  different  material. 

When  the  slit  cross  section  is  rectangular,  an  approximate 
solution  is  developed  in  Chapter  7  based  on  assumptions  concerning 
the  equivalent  magnetic  current.  This  approximate  solution  has  the 
advantage  of  yielding  a  greatly  simplified  numerical  solution,  accurate 
when  the  material  filling  the  slit  is  dense.  A  material  is  said  to  be 
dense  when  the  intrinsic  wavelength  in  the  material  is  much  smaller 
than  that  of  free  space. 
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Chapter  2 


PROBLEM  FORMULATION 

2.1.  Derivation  of  Operator  Equations 

The  problem  to  be  considered  is  shown  in  Fig.  1,  where 
electromagnetic  transmission  occurs  through  a  slit  of  arbitrary 
cross  section  in  a  perfectly  conducting  plane  of  finite  thick¬ 
ness.  The  slit  cross  section  is  composed  of  the  straight  lines 
T1  and  r3,  which  are  the  two  aperture  faces,  and  the  lines  and 
r4,  which  are  arbitrary  in  their  specification.  The  characteristic 
dimensions  of  the  slit  cross  section  are  shown  in  Fig.  1.  This 
problem  is  strictly  two-dimensional  in  that  everything  is  invar¬ 
iant  in  the  z  direction.  A  time  dependence  of  e^Wt  is  assumed 
throughout.. 

In  the  original  problem,  (E^+E3,  H*+Ha)  is  the  total  field 

b  b  c  c 

in  region  a,  ( E  ,  H  )  is  the  total  field  in  region  b,  and  (E  ,  H  ) 

is  the  total  field  in  region  c.  The  field  (JS*,  H*)  is  that  which 

would  exist  if  the  sources  (J*,  M*)  were  to  radiate  into  unbounded, 

homogeneous  space  with  (p  ,  e  )  everywhere.  The  equivalence  principle 

SL  Si 

[23,24]  is  used  to  isolate  the  three  regions  by  postulating  equiv¬ 
alent  sources  to  support  the  fields  in  these  regions. 

The  equivalent  sources  for  region  a  consist  of  an  electric 
current  sheet  £  which  exists  everywhere  on  the  plane  at  x  =  0  and  a 
magnetic  current  sheet  M^  which  exists  only  over  the  aperture  face 
This  current  distribution  radiates  into  unbounded  space  with  (ha,  £a) 
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Fig.  1.  Original  problem:  Slit  of  arbitrary  cross  section  cut  in 
a  conducting  plane  of  finite  thickness. 


everywhere  and  in  the  presence  of  the  sources  (J  ,  M1)  to  give  the 
fields  (E*  +  j£a,  +  Ha)  to  the  left  of  the  x  =  0  plane  and  zero 
field  to  the  right.  This  situation  is  shown  in  Fig.  2a.  The  con¬ 
dition  that  the  tangential  magnetic  field  is  zero  just  to  the  right 
of  the  plane  at  x  =  0  is  written  as 

l£a(Ja,  M1)  =  -  ]£  (2-1) 

The  operator  L^q  gives  the  magnetic  field  of  sources  radiating  in 
unbounded  space  with  (u^.  e  )  everywhere.  The  permittivity  eq  and 
permeability  u  are,  in  general,  complex  numbers.  The  index  q  stands 

q 

for  a,  b,  or  c  and  the  subscript  t  denotes  the  tangential  component 
found  by  the  usual  -n  x  n  x  operation.  This  magnetic  field  operator 
is  defined  explicitly  in  terms  of  electric  and  magnetic  current  sources 
by  the  equation  [24,  p.  130] 


n  x  Lhq(J,M)  =  ±  y  J(r)  -  L 


J 


n  x  J(r')  x  V  H^2) (kq | r-r ' | )dt ' 


1  \ 


M(r')  H(2)(k  I r-r 1 | )dt '  (2-2) 

- o  q - 


Here  the  domain  of  the  integrals,  which  is  usually  given  by  J  and 

4 

M,  is  restricted  to  the  contour  C  ■>  U^r  in  anticipation  of  the  fact 
that,  when  Eq.  (2-1)  is  enforced  at  F^,  not  all  of  the  electric  cur¬ 
rent  on  the  plane  x  =  0  contributes  to  the  tangential  magnetic  field 
there.  The  minus  sign  is  used  when  the  field  point,  £,  is  on  the  — n 
side  of  C(q=a  or  c)  and  the  plus  sign  is  used  when  q=b.  The  elemental 
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1 


arc  length  is  it',  1_  is  a  unit  dyad,  k  =  ,  r|  =  /y  /e  , 

(2) 

Hq  is  the  Hankel  function  of  second  kind,  order  zero,  and  _r '  de¬ 
termines  a  point  on  the  source  distribution. 

The  equivalent  sources  for  region  c  are  given  by  an  electric 
c 

current  sheet  £  which  exists  everywhere  on  the  plane  at  x  =  d  and  a 

3 

magnetic  current  sheet  M  which  exists  only  at  the  aperture  face 

This  current  distribution  radiates  into  unbounded  space  with  (y  ,  e  ) 

c  c 

everywhere  to  produce  the  field  (EC,  HC)  to  the  right  of  the  plane 
at  x  =  d  and  zero  field  to  the  left.  This  situation  is  shown  in  Fig.  2b. 
The  condition  that  the  tangential  magnetic  field  is  zero  just  to  the  left 
of  the  plane  at  x  =  d  is  written  as 

L^C(JC,  M3)  =  0  (2-3) 

t  ~  — 

he 

where  is  obtained  from  Eq.  (2-2). 

The  true  solution  to  the  problem  in  Fig.  1  requires  that  the 

tangential  electric  and  magnetic  fields  be  continuous  at  the  aperture 

faces  ri  and  Hence  the  equivalent  sources  for  region  b  consist  of 

1  3 

magnetic  current  sheets  -M  on  and  -M  on  along  with  an  electric 

b  _  U  b  cL 

current  -J  which  exists  on  the  contour  C  =  where  £  =  £  on 

be  * 

r  and  J_  =  J_  on  This  distribution  of  currents  on  the  contour  C 

radiates  into  unbounded  space  with  (h-^,  £^)  everywhere  to  produce  the 

field  (E^*3,  fl*5)  inside  C  and  zero  field  outside.  This  situation  is 

shown  in  Fig.  2c.  The  condition  that  the  tangential  magnetic  field  be 

zero  just  outside  the  contour  C  is  written  as 

n  x  Lhb(-Jb,  -M1  -M3)  =  0  (2-4) 
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where  Lbb  is  obtained  from  Eq.  (2-2).  Alternatively,  the  condition 
that  the  total  tangential  electric  field  be  zero  just  outside  the 
contour  C  is  written  as 

L^b(-Jb,  -M1  -M3)  =0  (2-5) 

6CJ 

where  the  operator  L  gives  the  electric  field  due  to  sources 
radiating  into  unbounded  space  with  (y  ,  e  )  everywhere  and  is 

q  q 

defined  by  the  equation  [24,  p.  130] 


n  x  Leq(J,M)  =  +  f  M(r)  + 


-  1 


nxM(r')  x  V  H^2)(kq|r-r' |)dt’ 


k  n 

-  -3J-  fix  (1  +  -±r  V  V  •) 
4  _  -  k2  -  - 

q 


J(r')  H^2)(kq|r-r’t)dt' 


(2-6) 


where  the  same  remarks  hold  as  those  following  Eq.  (2-2)  except  that 
the  plus  sign  is  used  when  q  =  a  or  c  and  the  minus  sign  is  used  when 
q  =  b. 

An  alternative  electric  field  equation  could  also  have  been 
written  in  place  of  Eqs.  (2-1)  or  (2-3).  This  is  not  done,  however, 

SL  C 

because  the  electric  currents  J  and  J_  over  the  infinite  portions  of 
the  planes  at  x  =  0  and  x  =  d  contribute  to  the  tangential  electric 
field  on  and  respectively.  These  electric  currents  do  not  con¬ 
tribute  to  the  tangential  magnetic  field  on  and  F^,  and  hence  do 
not  appear  in  Eqs.  (2-1)  and  (2-3)  when  they  are  enforced  on  and 

r3  respectively.  The  necessary  unknown  equivalent  sources,  then, 
13  b 

consist  of  M  ,  M  ,  and  J[  which  occupy  a  finite  domain  given  by  the 
closed  two-dimensional  contour  C. 

10 

. — ■  * —  - - ----- 


The  solution  to  the  operator  equations  (2-1),  (2-3),  and 
either  (2-4)  or  (2-5)  will  be  unique  if  the  sets  of  homogeneous 
equations  consisting  of 

L^a(Ja,  M1)  =  0  (2-7) 

and  either  (2-3)  and  (2-4)  or  (2-3)  and  (2-5)  have  only  the  trivial 
solution.  Thus  it  remains  to  show  that  if  the  sources  are  removed 
from  the  problem  in  Fig.  1,  the  fields  everywhere  will  collapse  to 
zero. 

Equation  (2-7)  states  that  the  tangential  magnetic  field 
equals  zero  just  to  the  right  of  the  plane  at  x  =  0  due  to  J_  and  M 
radiating  in  unbounded  space  with  (y  ,  £  )  everywhere.  This  current 

3  3 

SL  cL 

distribution  also  creates  the  fields  (E  ,  H  )  everywhere  to  the  left 

of  the  plane  at  x  =  0.  This  situation  is  shown  in  Fig.  3a.  Equations 

(2-4)  and  (2-5)  state  that  the  tangential  electric  and  magnetic  fields 

1  3  b 

are  zero  just  outside  C  due  to  -M  ,  -M  ,  and  -J  radiating  in  unbounded 

space  with  (y,  ,  £  )  everywhere.  This  implies  that  the  electromagnetic 

field  everywhere  outside  C  is  zero  while  the  field  inside  C  is  given  by 

(E^,  H^).  This  situation  is  shown  in  Fig.  3b  and  it  is  noted  that 

either  of  Eqs.  (2-4)  or  (2-5)  is  sufficient  to  produce  this  situation. 

Equation  (2-3)  states  that  the  tangential  magnetic  field  equals  zero  just 

c  3 

to  the  left  of  the  plane  at  x  =  d  due  to  J_  and  M  radiating  in  unbounded 

space  (y^,  ec)  everywhere.  This  current  distribution  also  creates  the 
c  c 

fields  (E  ,  H  )  everywhere  to  the  right  of  the  plane  at  x  =  d.  This 
situation  is  shown  in  Fig.  3c. 


zero 

field 


zero 

field 


(V£a) 


b,  Hb 


(Vv,e.  ) 


—  J^^^,plane  x  =  0 

Fig.  3a.  Situation  of  Eq.(2-7) 
Ja  =  Jb  at  r_ . 


-J 

contour  C 


Fig.  3b.  Situation  of  Eqs.(2-4) 
and  (2-5),  -Jb  =  -Ja  on 
F1  and  ~J  =  -J0  on  r 
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Fig.  3c.  Situation  of  Eq.(2-3) 

JC 

=  J  at  r„. 
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conductor 
\  \  \  \  \  \ 
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(ub’eb} 


perfect  \ 
\  conductor 


Fig.  3d.  Composite  situation 


Now  the  composite  situation  shown  in  Fig.  3d  is  considered 

where  the  fields  of  Figs.  3a  -  3c  have  been  superposed.  From  Fig.  3a, 

the  tangential  component  of  _E  is  zero  everywhere  at  x  =  0  except  on 

c 

r  and  from  Fig.  3c,  the  tangential  component  of  is  zero  every¬ 
where  at  x  =  d  except  on  From  Fig.  3b,  the  tangential  component 

of  is  zero  at  the  perfectly  conducting  contours  Y^  and  T , .  Figures 
3a  and  3b  indicate  that  there  are  no  equivalent  sources  at  in  Fig.  3d 
and  Figures  3b  and  3c  indicate  that  there  are  no  equivalent  sources  at 
in  Fig.  3d.  The  conclusion  is  that  the  tangential  electric  field 
is  zero  everywhere  on  the  boundaries  of  the  perfect  conductors  in  the 
composite  situation  of  Fig.  3d  and  also  that  the  region  outside  the 
perfect  conductors  is  source  free.  Thus  the  fields  in  Figs.  3a-3c  are 
collapsed  to  the  null  field  and  hence  the  solution  to  Eqs.  (2-1),  (2-3), 
and  either  (2-4)  or  (2-5)  is  unique. 

2,3.  Specification  of  Slit  Cross  Section,  Basis  Functions,  and 
Symmetric  Product 

The  contour  C  which  determines  the  slit  cross  section  is 
approximated  by  a  number  of  straight  line  segments  each  of  length  AC 

n 

for  integers  n=l ,2 , . . . ,N^-1.  This  is  shown  in  Fig.  4.  The  integers 

N^,  i=l,2,3,4,  or  5  are  assigned  to  the  beginning  and  end  points  of 

each  T.  with  N,  =  1  so  that  each  T.  is  broken  up  into  N.,,  -  N.  seg- 

i  1  l  l+l  i 

ments.  Each  segment  AC  has  a  transverse  directed  normal  n  and 

n  — n 

tangent  £  such  that 
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Fig.  5.  Typical  segments  AC  and  AC  .  R  is  defined  as  a  vector 

m  n  — m  ,n 

from  the  midpoint  of  AC  to  the  midpoint  of  AC  . 


z 


(2-8) 


£  x  ft  = 

— n  — n 

It  is  useful  to  introduce  the  parameter  t  which  represents  the  arc 
length  along  C  from  the  origin  to  any  point  on  C.  Thus  the  sub¬ 
scripted  values  of  t  are  given  by 

n-1 

t  =  l  AC  (2-9) 

i=l 

for  n=2,3,...,  and  t^  =  0.  Two  typical  straight  line  segments  are 

shown  in  Fig.  5  where  the  vector  R  is  also  defined. 

— m,n 

For  each  polarization,  both  transverse  (£)  and  longitudinal 
(z)  directed  current  sources  are  present.  The  transverse  currents 
have  a  charge  associated  with  them  and  it  is  desirable,  though  not 
necessary,  that  their  representation  be  differentiable.  The  z  directed 
currents  have  no  such  requirements  but  sometimes  become  unbounded  at 
edges.  The  longitudinal  expansion  functions,  then,  should  not  be  con¬ 
strained  to  be  continuous  at  edges.  For  the  sake  of  simplicity,  the 
edge  condition  [27,28]  will  not  be  built  in  a-priori  into  the  repre¬ 
sentation  of  the  unknown  currents.  With  these  considerations  in  mind, 
then,  two  different  sets  of  expansion  functions  are  defined  as 


and 


f  (t) 
— m 
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*m(t)  = 


(2-11) 


z 


for  t  <  t  <  t  ,, 
m  —  —  m+l 


for  t  elsewhere 


for  m=l,2,...,  Nc-1  and  with  t  =  t1T  .  . 

o  o  N^-l 

Lastly,  in  order  to  carry  out  the  method  of  moments  procedure, 
it  is  necessary  to  define  a  symmetric  product  as 


<A,  B> 


j  A  •  B  dt 
C 


(2-12) 


where  A  and  ]5  are  vector  functions  defined  on  the  two  dimensional  con¬ 
tour  C. 


2.4.  Reduction  of  Operator  Equations  to  Matrix  Equations 


The  method  of  moments  [25]  is  applied  somewhat  mechanically 
to  solve  either  one  of  two  possible  sets  of  simultaneous  operator  equa¬ 
tions.  These  consist  of  the  sets  (2-1),  (2-3),  and  (2-4)  or  (2-1),  (2-3), 
and  (2-5).  The  two  different  polarizations  are  considered  separately 
in  Chapters  3  and  4.  The  first  step  is  to  assume  that  the  unknown  cur¬ 
rents  can  be  represented  by  a  linear  combination  of  vector  basis  func¬ 
tions  defined  on  the  two-dimensional  contour  C.  These  may  be  chosen  from 
the  sets  defined  in  Section  2,3.  Thus  the  equivalent  currents  are 
represented  in  general  by 


- 1 

v1 

n 

%n(t) 

(2-13) 

n 

- 1 

v3 

n 

-3n<C> 

(2-14) 

n 

=  I 

n 

I 

n 

^(t) 

(2-15) 
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where  V  ,  V  ,  and  I  are  unknown  complex  coefficients  to  be  determined, 
n  n  n 

The  functions  a^,  $  are  vector  basis  functions  defined  on  1^, 

and  C  respectively  and  will  be  chosen  from  the  sets  (2-10)  and  (2-11) 
depending  on  the  polarization  considered. 

The  next  step  is  to  substitute  Eqs.  (2-13)  -  (2-15)  into  the 
operator  equations  and  "test"  each  equation  on  that  portion  of  the  con¬ 
tour  C  over  which  it  is  valid.  This  is  done  by  taking  the  symmetric 
product  of  the  operator  equation  with  "testing"  functions  defined  on  the 
contour  C.  It  is  desirable  to  choose  testing  functions  which  are  in  the 
range  of  the  operator.  The  operator  equations  have  been  written  in  such 
a  way  as  to  allow  the  use  of  the  basis  functions  for  the  testing  procedure. 
This  is  computationally  convenient  in  that  several  matrices  become  sym¬ 
metric.  Thus  Eq.  (2-1)  is  tested  with  a,  at  T,  ,  Eq.  (2-3)  is  tested  with 

— lm  1 

a.  at  T. ,  and  Eqs.  (2-4)  and  (2-5)  are  tested  all  along  C  with  0  .  This 
— Jm  j  °  — m 

procedure  results  in  two  possible  sets  of  simultaneous  equations  which  are 
written  in  matrix  form  as 

[Y3]  [  0  1 

[  0  ]  [YC] 

[Yhl]  [Yh3] 
or 

[Y3]  [  0  ] 

[0]  [YC] 

[Y61]  [Ye3] 

The  matrix  [  0  ]  stands  for  a  matrix  of  zeros  and  5  stands  for  a  column 
vector  of  zeros.  As  shown  earlier,  either  (2-16)  or  (2-17)  provide  a 


[U1] 

I1 

[U3] 

= 

3 

(2-16) 

[Th] 

3_ 

[U1] 

"v1" 

¥ 

tu3] 

v3 

= 

3 

(2-17) 

[T6] 

V 

3 
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unique  solution  to  the  original  problem.  The  matrix  elements  are  given 
by  the  general  formulas: 


Ya  =  -  k  n  <a,  ,  Lha(0,  a,  )> 
mn  o  o  — lm  — t  -In 


(2-18) 


mn 


mn 


-  Vo  *3n>> 

k 

- 7T  <a  ,  n  *  8  > 

2  — lm  —  — n 


(2-19) 


(2-20) 


mn 


V  <&-  ,  fi  x  6  > 

2  3  m  —  — n 


(2-21) 


=  -  k  <2  >  fi  X  Lhb(0.  fl1n)> 

mn  o  o  — m  —  —  -  In 


(2-22) 


Yb3  =  -  k  x  Lhb(0,  a.  )> 

mn  o  o  — m  —  —  — ’  — 3n 


(2-23) 


=  -  k  <0  ,  fi  x  Lhb(e  ,  0)> 

mn  o  ~m  —  —  ^-n  — 


(2-24) 


'f61  -  -  <4-  if  <0,  aln)> 


mn 


(2-25) 


Ye3  =  -  k  <8  ,  Leb(0,  a_  )> 
mn  o  — m  — t  —  — 3n 


(2-26) 


T6  =  -  —  <8  ,  Leb(g  ,  0)> 
mn  r)  -m  — t  — n’  — 
o 


(2-27) 


The  linearity  of  the  operators  Leq  and  Lhq  has  been  used  in  the  above. 
Scaling  has  also  been  done  with  k^,  the  wavenumber  of  free  space,  and 
0  »  the  impedance  of  free  space,  so  that  the  computed  matrix  quantities 
are  dimensionless.  The  minus  signs  are  kept  in  (2-18)-(2-27)  so  that 
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Fig.  6.  Network  coupling  interpretation  of  partitioned  matrix  equations. 


Chapter  3 

TRANSVERSE  ELECTRIC  (TE)  CASE 


3.1.  Expansion  of  the  Equivalent  Currents 


For  the  TE  polarization,  there  is  only  a  z  component  of  magnetic 
field  and  transverse  components  of  electric  field.  The  equivalent  mag¬ 
netic  currents  are  z  directed  and  satisfy  no  boundary  conditions  on  r 
and  T  .  The  equivalent  electric  currents  are  transverse  directed  and  are 
continuous  all  along  the  contour  C.  With  this  in  mind,  the  functions  of 
Eq.  (2-11)  are  chosen  for  and  of  (2-13)  and  (2-14)  and  the  func¬ 
tions  of  Eq.  (2-10)  are  chosen  for  £  of  (2-15).  The  unknown  currents  are 


then  expanded  as 


and 


1  c  1 

M  =  l  Vn  JB^t) 


VI 

m3=  l  v3^(  t) 

n=N^ 


J  =  I  I  f  (t) 

n=l  n  ^ 


(3-1) 

(3-2) 

(3-3) 


The  formulas  for  the  matrix  elements  of  (2-16)  and  (2-17)  may  now  be 
obtained  from  Eqs.  (2-18)  -  (2-28).  They  are  given  as  follows  where 
the  ranges  of  the  indices  m  and  n  are  specified: 


1  £  m  <  N  -1 

kono  "V  .  <  <  . 

1  £  n  <_  N^-l 


(3-4) 


YC  =  -  k  n 
mn  o  o 


N,  <  m  <  N  -1 
3  —  —  4 

N,  <  n  <  N  -1 
3  —  —  4 


(3-5) 
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mn 


n  x  f  > 
—  — n 


mn 


n  x  f  > 
—  — n 


rhl 


hb , 


-  k0n0  <4.  A  *  I  C°,  s„>> 


f"3 

mn 


hb. 


=  -  k  n  <f  ,  n  x  L  (0,  g  ) 
o  o  -m  —  -  —  °txj 


mn 


k 

o 


n  x  Lhb(f  ,  0)> 
— n  — 


fel 

mn 


eb , 


k  <f  ,  L7  (0 ,  g  )> 
o  -m  — t  —  ’ 


>re3 

mn 


eb  . 


*„  <4.  4  <0,  s„)> 


mn 


-  —  <f  , 
n  — m 

o 


0)> 


1  £  m  <  N2-l 
1  <  n  £  N5-l 

N,  <  m  <  N.-l 
1  £  n  £  N  -1 

1  £  m  £  N  -1 
1  £  n  £  N  -1 

1  £  m  £  N  -1 

N_  <  n  <  N.-l 

J  —  —  4 

1  £  m  £  N  -1 

1  £  n  £  N  -1 

1  £  m  £  N5~l 
1  £  n  £  N  -1 

1  £  m  £  N  -1 

N.  <  n  <  N  -1 

i  ~  —  4 

1  £  m  £  N  -1 

1  £  n  £  N  -1 


The  non-zero  elements  of  the  excitation  vector  are  given  by 


(3-6) 

(3-7) 

(3-8) 

(3-9) 

(3-10) 

(3-11) 

(3-12) 

(3-13) 


m 


kono  "V 


1  £  m  £  N2-l  (3-14) 
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3.2.  Formulas  for  TE  Matrix  Elements 


Equations  (3-4)  -  (3-13)  define  the  elements  of  the  matrices  in 
Eqs.  (2-16)  and  (2-17).  Explicit  formulas  will  be  derived  here  for 
these  symmetric  products  using  the  definitions  of  the  operators  given 
by  Eqs.  (2-2)  and  (2-6)  and  the  expansion  functions  (2-10)  and  (2-11). 
Equations  (3-4)  and  (3-5)  may  be  computed  by  examining  the  inner  pro¬ 
duct 


ha . 


P1  ■  -  Vo  <V  it  S„» 


k  n  k 
o  o  a 

4n 


"m+l 


n+1 


^(t)*  fix  fix  I  Sn(t,)H^2)(ka|rm(t)  -r^(t)|)dt,dt  (3-15) 


m  n 

where  r  (t)  denotes  a  vector  from  the  origin  to  a  point  on  AC  for 
— m  m 

2  2 

m=l,2 . Hg-1.  Using  the  variables  u  =  t  ,  u’  =  t'  ,  and 

m  n 

R  defined  in  Fig.  5,  one  may  write  Eq.  (3-15)  as 
— m,  n 

kono  Ym  Yn  f  f  (2).iYm  Yn  ip  +  v  R  |  x  H  'a 

H  (  hr  ut r-  u '  t  +  k  R  )du  du 

o  2  — m  2  — n  a— m ,  n 1 


4  k  r|  4 
a  a 


J 

-1  -1 


P1  = 


if  m  4  n 


(3-16) 


k  n  Ym 

o  o  m 


4  k  n  2 
a  a 


[a  (u+1))  +a  (-f  (l-u))]du 


if  m  =  n 


where  Y  =  k  AC  ,  Y  =  k  AC  ,  and  the  function  a  is  defined  by 
m  a  m  n  a  n 


a(z) 


z 

(2) 

IT  ;  (v)dv 
o 

0 


(3-17) 


This  function  may  be  computed  using  Struve  functions  [29-11.1.7].  The 
integrals  in  (3-16)  may  be  approximated  readily  by  Gaussian  quadrature 
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m-1  'n 
16 


J  J  (2  +  Ho2)<|JTi  u  4-1  -  f  "'4  +  kA-i,nl>d“'i» 


if  m-1  ^  n 

(3-21) 


(f +  f)[a  (y  (u+1))  +  a(y  (l-u))]du 


if  m-1  = 


^Y„  Y.  1  1 


m  n 


/ZSL  +  i\  «(2)/i  m  ^  f 

16  J  I  (  ?  +  Y)  H„  (|—  u  1_  "  —  u 


-1  -1 


T+2>  Ho  (lTu4-tu'4  +  kb\,nDdu,du 


S2  = 


Ym  r 


if  m  ^  n 

(3-22) 


(y  +  j)[a  (y  (u+1))  +  a  (y  (l-u))]du 


-1 


if  in  =  n 


For  Eq.  (3-10),  the  magnetic  field  due  to  a  transverse  directed 
electric  current  is  needed  and  using  (2-2)  with  M  =  0,  one  obtains 


P,  = 


hb. 


k  <f  ,  n  x  L  (f  ,  0)> 
o  -m  —  —  '-n’  — 


k  f1 

o 


f  (t)  •  f  (t)dt 


m-1 


k  f1 

o  1 


n+1 

r 


4(t)  *  nxf^  (t>)  x  V  H^2)(kb|r-r,|)dt’dt  (3-23) 


m-1 


n-1 


where  r  is  a  vector  from  the  origin  to  a  point  on  AC  ,  U  AC  and  r' 

m-1  m  — 


is  a  vector  from  the  origin  to  a  point  on  AC  u  AC  .  The  above 


integral  is  broken  up  into  four  parts  and  the  symmetric  product  is 


written  as 


=  -j—  [S^(m-l,n-l,l,l  )  +  S^(m-l,n,l,-l) 


+  S^(m,n-l,-l,l)  +  S^(m,n,-1,-1) ] 


(3-24) 


The  function  is  defined  by 


(m,n,p,q) 


1  1 


-1  -1 


(P  Y  +  |)  (Q  |  +  |)  4>(m,n)du'du  if  m^n 


(3-25) 


m  ,1  ,  pq> 

"  —  (2  +  6  )  lf  m=n 


where 


n  •  ft  ,9. 

<t>  (m,n)  =  — - =^H(2)(|6  I) 

|r  I  1  _m»n 

— m,n 


(3-26) 


ft  =  k,  R  +  —  u£ - —  u'£ 

— m,n  trTn,n  2  — m  2  — n 


(3-27) 


A  similar  transformation  of  variables  as  was  done  for  Eq.  (3-19)  has 
been  used  as  well  as  the  identity 


V  H^2)(kb|r-r' |)  =  -^ - -  Hj2)(kb|r-r'|) 

I r-r ' I 


(3-28) 


Equations  (3-11)  and  (3-12)  are  computed  using  the  electric 


field  operator  defined  in  Eq.  (2-6),  from  which  one  obtains 
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p5  =  -  ko  <4’  V<°’  Sj> 


<Km,n) 


A  A 

n  *  R 

~~m  — m,n  u(2) 


(3-33) 


where  R  is  defined  in  Eq.  (3-27). 

— m,n 

Finally,  to  obtain  Eq.  (3-13),  one  must  compute  the  following 
symmetric  product: 


P,  =  -  —  <f  ,  LSb(f  ,  0)> 

6  n  ~m  — t  -n  — 

O 

k  k  H  m+1  rn+1 

=  "  ~T„—  f  (t)  •  nxfix  (l  +  i  V  V  •)  f  (t')HU;(k  |r-r'l)dt'dt 

°  i  .  N.  i  , 

"-1  “-1  (3-34) 


where  r  is  a  vector  from  the  origin  to  a  point  on  AC  ,  u  AC  and  r '  is  a 
—  m-I  m  — 

vector  from  the  origin  to  a  point  on  AC  U  AC  .  The  technique  of  inte- 

n-1  n 

gration  by  parts  is  applied  twice  to  the  term  containing  the  V  V  •  operation 
and  the  result  is  given  by: 


kokbnb  P*1  fn+1 


'  6  4n 


df  (t)  df  (f) 

!4(t)  ■io(t'>  -  :r~dt  >(kbl-"i'l)dt'dt 


t  ,  t  1 
m-1  n-1 


(3-35) 


This  integral  is  broken  up  into  four  parts  and  after  similar  trans¬ 
formations  as  used  in  (3-19),  one  obtains 


P6  =  kf'rT  n-1, 1,1)  +  S6(m-1,  n,  1,  -1) 


+  (m,  n-1,  -1,  1)  +  S6(m,  n,  -1,  -1)} 


(3-36) 


where  the  function  S,  is  defined  by 

D 
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1  1 


uH<2)(u)du  =  z  hJ2)(z)  -  ZL  (3-38) 

0 

3.3.  Formulas  for  TE  Excitation  Vector 

The  elements  of  the  non-zero  excitation  vector  are  given  by  Eq.  (3-14). 
The  incident  magnetic  field  is  z  directed  and  is  adjusted  so  that  the  inci¬ 
dent  electric  field  is  equal  to  unity  and  has  zero  phase  at  the  center  of 
aperture  face  with  respect  to  the  coordinate  system  in  Fig.  4.  Plane 
wave  and  line  source  excitations  are  considered  separately.  For  plane 
wave  excitation,  the  incident  magnetic  field  is  given  by 
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Chapter  4 


TRANSVERSE  MAGNETIC  (TM)  CASE 


4.1.  Expansion  of  the  Equivalent  Currents 

For  the  TM  polarization,  there  is  only  a  z  component  of  electric 
field  and  transverse  components  of  magnetic  field.  The  equivalent  mag¬ 
netic  currents  are  transverse  directed  and  must  satisfy  the  boundary 
condition  that  they  are  zero  at  the  ends  of  and  F^.  The  equivalent 
electric  currents  are  ^d^re^ted  an4  satisfy^o-  botmiary  cendi^ons  *bn 
the  contour  C.  In  fact,  it  is  expected  that  becomes  unbounded  at 
the  ends  of  F^  and  when  region  b  is  lossless  since  it  is  the  actual 
electric  current  induced  on  these  surfaces  in  the  original  problem. 

With  this  in  mind,  the  functions  of  Eq .  (2-10)  are  chosen  for  the  and 

of  (2-13)  and  (2-14)  while  the  functions  of  Eq.  (2-11)  are  chosen  for 
of  (2-15) .  The  unknown  currents  are  thus  expanded  as 


M1  =  l  v]  f  (t)  (4-1) 

n=2  n 

V1 

M 3  =  ?  V3  f  (t)  (4-2) 

—  “  ,  i  n  — n 

n=N^+l 
b  >'5"1 

i  -  \  W1’  <4-3> 

n=l 


The  formulas  for  the  matrix  elements  of  Eqs.  (2-16)  and  (2-17)  are 
obtainable  from  Eqs.  (2-18)-(2-28) .  They  are  given  as  follows  where 
the  ranges  of  the  indices  m  and  n  are  now  specified: 
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The  non-zero  elements  of  the  excitation  vector  are  given  by 


El'  -  * 


I1  =  k  n  <f  ,  Hx> 
m  oo  — m  — t 


2  <  m  <  N2-l 


(4-14) 


4.2.  Formulas  for  TM  Matrix  Elements 


Equations  (4-4)  -  (4-13)  define  the  elements  of  the  matrices 
in  Eqs .  (2-16)  and  (2-17).  Explicit  formulas  will  be  derived  here 
for  these  symmetric  products  using  the  definitions  of  the  operators 
given  by  Eqs.  (2-2)  and  (2-6)  and  the  expansion  functions  (2-10)  and 
(2-11). 

For  Eqs.  (4-4)  and  (4-5)  it  is  necessary  to  know  the  magnetic 
field  due  to  a  transverse  directed  magnetic  current.  This  is  obtained 
by  duality  from  the  electric  field  due  to  a  transverse  electric  current 
which  was  used  in  Eq.  (3-13)  for  the  TE  case.  Thus  one  may  write 


k  n 

(S6(m-1,  n-1 ,  1,  1)  +  S6(m-1,  n,  1,  -1) 


+  S6(m,  n-1,  -1,  1)  +  S6(m,  n,  -1,  -1)}  (4-15) 


where  the  function  S,  is  defined  in  Eq.  (3-37)  and  k  replaces  k,  . 

b  a  b 

Equations  (4-6)  and  (4-7)  are  straightforward  and  one  may 
easily  write 
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- - - - 


q2  =  -  t  <4*  £  *  V 


k  AC 
o  n 


0 


for  m-1  =  n  or 


for  m  otherwise 


m  =  n 


(4-16) 


Equations  (4-8)  and  (4-9)  require  rhe  magnetic  field  due  to  a 
transverse  directed  magnetic  current.  Thus  the  inner  product  is 
written  as 

«3  ‘  -  Vo  <V  a  *  ihb«.  V* 


kbkonQ  ^ 

4nb  J 


n+1 


”n-l 


fn(f)H^2)(kb|r-r’|)dt'dt 


kono 

4kb\ 


rm+1 


/o\ 

V-j  f^t^H^'OcJr-r’ |)dt’dt  (4-17) 


n-1 


where  £  is  a  vector  from  tne  origin  to  a  point  on  AC  and  r'  is  a  vector 

m  — 

from  the  origin  to  a  point  on  AC^  U  AC^.  After  some  manipulation,  one 
may  write  the  above  in  the  form 


k  n  6 
j  =  -2-2.  y  s 

3  kbnb  i-1  1 


•  (4-18) 


where  each  is  defined  by 
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(4-23) 
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(4-24) 


4^  °<V 

n 


In  the  above,  Y  =  k.  AC  . 

m  b  m 


if  m  =  n  or  n+1 


For  Eq.  (4-10),  the  magnetic  field  due  to  a  z  directed  electric 
current  is  necessary  and  one  obtains 
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m  Cn 

where  r_  and  r'  are  vectors  from  the 


origin  to  points  on  AC  and  AC 


respectively.  The  above  is  written  in  the  form 

r  K 


Q4  = 


AC 

£  m 


ko  1  Vo 
tb 4J  4 


if 


m  =  n 


-1 


1  1  - 

r  n  •  R 

~n  ~m,n  (2') 
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(4-26) 


where  Ym  ^A^  and  is  defined  by  Eq.  (3-27). 


For  Eqs .  (4-11)  and  (4-12)  it  is 


field  due  to  a  transverse 


necessary  to  know  the  electric 


magnetic  current.  Thus 
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where  r  is  a  vector  from  the  origin  to  a  point  on  AC  and  r'  is  a  vector 

m  — 

from  the  origin  to  a  point  on  AC^_^  U  AC^.  This  equation  is  very  similar 
to  Eq.  (3-29)  and  after  some  manipulation,  one  obtains 

k 


Q5  -  ^  (S,  +  s2) 


(4-28) 


where  and  S2  are  defined  by 


if  m  =  n-1 


{~2  +  2”)  <Km,n-l)du'du  if  m  ^  n-1 


(4-29) 


if  m  =  n 


(“  -y"  +  -j)  <f(m,n)du'du  if  m  ^  n 


-1  -1 


(4-30) 


and  (^(nijn)  is  defined  in  Eq.  (3-26).  In  the  above,  y  =  k.AC  . 

m  Dm 

Lastly,  for  Eq.  (4-13),  it  is  necessary  to  have  the  electric  field 
due  to  a  z  directed  electric  current.  This  is  found  simply  from  the  mag¬ 
netic  field  due  to  a  z  directed  magnetic  current  used  in  Eq.  (3-15).  Hence 
one  easily  obtains 


Q6  =  ‘F 


,  eb 


%’  h  (v 


Tpi 


(4-31) 


where  P  is  defined  by  Eq.  (3-16) 
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4.3.  Formulas  for  TM  Excitation  Vector 


The  elements  of  the  non-zero  excitation  vector  are  given  by 
Eq.  (4-14).  The  incident  electric  field  is  z  directed  and  is  adjusted 
so  that  it  is  equal  to  unity  and  has  zero  phase  at  the  center  of  aper¬ 
ture  face  r  with  respect  to  the  coordinate  system  in  Fig.  4.  Plane 
wave  and  line  source  excitations  are  considered  separately. 

For  plane  wave  excitation,  the  incident  magnetic  field  has  a 
component  tangential  to  given  by 


i  _  cos  j)1  e"jka[x  COS  4,1  -  (y  "  T5  Sin  ^ 

y  "  n„ 


where  the  angle  of  incidence  4>  is  measured  from  the  negative  x  axis. 
Substituting  this  into  Eq.  (4-18),  one  obtains 


w  i 

,  „  .  a.  .  ,i  -jk  AC  ,sin  <p 

■  k  1  a  1  jk  (t - x-)sin  4>  a  m-1 

tol  o  o  cos  ®  J  a  m  2  re  -1 

X  ~  ■;  .  e  I - - - 

m  k  n  .  .  .1  „  .  T 

a  a  j  sm  q>  ik  AC  ,sin  d> 

a  m-1 


(4-32) 


jk  AC  sin  61 
a  m 


jk  AC  sin  cj)1 
am 


(4-33) 


If  an  electric  line  source  of  strength  I1  amps  is  placed  at  the 
coordinates  (xg,  y  )  in  region  a,  then  the  tangential  component  of  mag¬ 
netic  field  at  F^  is  given  by 


H 


i 

y 


x 

s 


*s  +  (y  -  ys)2) 


+  (y  -  ys) 


2 


-  y.)  > 


where  I  has  been  adjusted  to 


(4-34) 


Wo  '<V  xs  +  <f>  > 
Substituting  Eq.  (4-34)  into  (4-14),  one  obtains 
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Chapter  5 

MEASURED  QUANTITIES 


Once  the  magnetic  currents  on  the  aperture  faces  are  determined, 
several  transmission  characteristics  of  the  slit  may  then  be  computed. 
In  each  following  section,  the  TE  case  is  treated  first,  followed  by 
the  TM  case. 


5.1.  Measurement  Vector 

For  the  TE  case  the  magnetic  current  sheet  at  is  z  directed, 

c 

hence,  to  measure  the  magnetic  field  H  at  a  point  in  region  c  as  shown 
in  Fig.  7,  a  magnetic  line  source  K  z_  is  placed  at  the  measurement  point. 
From  reciprocity  one  obtains 


H  •  K  z_  =  H  •  M  dy 


where  H  is  the  magnetic  field  due  to  K  z_  radiating  in  the  presence  of 

c  3 

a  complete  conductor  and  _H  is  the  magnetic  field  of  currents  M  radi¬ 
ating  in  the  presence  of  a  complete  conductor.  After  using  Eq .  (3-2) 

3 

for  M  ,  Eq.  (5-1)  becomes 


(5-1) 


■  l  K  ^  ■  t 


r  "-'<3 

3  J 


A  V1 

c  v 


j  Vn  F  8„(y)H‘2)(Kcly2-rml)dy  (5-2) 
n=N„  ■'w 


Here  the  coordinate  origin  is  now  in  the  center  of  and  the  above 
equation  may  be  written  in  matrix  form  as 


H 


c 

z 


-jk  r 
cm 


-*e  ~*'3 
I  V 


(5-3) 


where  the  tilde  (~)  denotes  transpose  and  I  is  a  TE  measurement  vector 
whose  elements  are  defined  as 


yirk  r  jk  r 
cm  cm 

IT' 


8„<y)Ho2)(kc|y2  -  t,')dy 


(5-4) 


for  n  =  N„ ,  N_+l,...,  N,-l.  If  r  »  X  so  that  far-field  measurements 
3  3  4  me 

are  of  interest,  the  large  argument  approximation  to  the  Hankel  function 
is  used  in  Eq.  (5-4)  to  obtain 


,ef 


w  /2 
c 


-w  /2 
c 


jk  y  sin  <J> 
gn(y)e  c  dy 


AC 


k  AC 


=  -  AC  e 
n 


jk  (h  +  — -i-)sin<J>  sin  ( — - —  sin  <)>) 
c  n  2.  i 


k  AC 

(  C2  -  sin  4>) 


(5-5) 


where  h  =  y-  -  \  A^  for  n  =  N  ,  N  +1 , . . . .N^-l .  This  can  be  con- 

n  i=N_  1 


veniently  written  in  terms  of  the  excitation  vector  of  Eq.  (3-40). 

For  the  TM  case,  the  magnetic  current  sheet  at  is  y  directed, 
hence  the  electric  field  E,C  in  region  c  is  measured  by  placing  an 

£/n 

electric  line  source  I  z  at  the  measurement  point  as  shown  in  Fig.  8. 
Again  from  reciprocity,  one  obtains 


E 


c 


Tl~ 

I  z 


dy 


(5-6) 
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Equation  (4-2)  is  substituted  for  M  and  the  above  is  rewritten  as 


ECI*  = 
z 


V1 


P  n=N.+l 
3 


y  v  f  (t)  •  h  dt 

u  n— n  — 


where  use  has  been  made  of  the  fact  that  f  (t)  on  T.  is  -y  directed. 

— T3  A 


Using  the  y  component  of  H  in  the  above,  one  obtains 


,  ,£  N  -1 
„  0  k  I  4  , 

,cT£  c  v  ,.3 

-I  = - 7TT-  X  )  V 

z  2j  m  J.1  n  J 


n=N  +1 

3  w 


f  H<2)(kJy£-rJ>dy 

fn(y)  - 


2  2 
*  +  (y  -  y  ) 

m  m 


where  r  =  x  x  +  y  y.  An  expression  for  E  may  now  be  written  in 
— m  m~  z 

matrix  form  as 


e‘  -  p-  f"  V3 

z  \  27rr 


m 


where  I  is  the  TM  measurement  vector  whose  elements  are  given  by 


x  /irk  r  jk  r 
_m  m  /  c  m  cm 

"  =  "  j  J  2j  6 


f  H^2)(k  |yy  -  r  |)dy 

2  f  (y)  1  c  - 


j  n  r 

W  /  X 

__c  \J 


2  2 

+  (y  -  y_) 

m  m 


for  n  =  N_+l,  N_+2,...,  N.-l.  If  r  >>  X  for  far  field  measurements 
3  3  4  me 


Eq.  (5-10)  becomes 


w 


w 


_mf 

I  =  -  cos i 
n 


£„(y)  e 


jk  y  sin  $ 


dy 


v. 

r 

k:  -  I 

». 


=  COS  <f) 


AC 

jk  (h - ^)sin  <t> 

2AC  e  c  n  2  k  AC 

- 2 - X -  sin  (  c9  n-  sin  $) 

j(kcACn  sin  <p) 


cos 


2AC  ,  e 
n-1 


AC 

jk  (h  +  — ?p4sin 
c  n  l 


j  (k  AC  sin  d>)  ‘ 
J  c  n-1  Y 


k  AC 

sin  (  -C  —  sin  0) 


(5-11) 


where  h  =  -  £  AC,  1  for  n=N  +1,  N  +2,...,N  -1.  If  each 


i=N3+l 

AC  =  AC  on  T  ,  then  (5-11)  becomes 
n  j 


L 


,mf 


cos  d>  AC  e 


jk  h  sin 
c  n 


k  AC 

sin  (— ^ —  sin  <J>) 


k  AC  sin  4> 
(“^ - ) 


5.2.  Transmission  Coefficient 


The  transmission  coefficient  of  the  slit  is  defined  by  the 


ratio 


T  = 


t3 

iN 


(5-12) 


(5-13) 


where  Pf3  is  the  time  average  power  transmitted  through  the  aperture  face 
r3  and  P^N  is  the  time  average  power  intercepted  by  the  aperture  face  P^ 
when  the  source  is  normally  incident.  If  the  time  average  power  trans¬ 
mitted  through  aperture  face  is  denoted  by  P  ,  then  when  the  slit 
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r 


is  filled  with  a  lossless  medium,  must  equal  P^*  This  may  be  used 
as  a  check  on  the  numerical  solution  when  the  imaginary  parts  of  the  com¬ 
plex  power  through  and  r  are  not  significantly  different. 

Consider  the  TE  case  first.  A  magnetic  line  source  is  placed 
w 

at  (xg,  -^)  in  region  a  to  produce  an  incident  field.  The  slit  inter¬ 
cepts  a  fraction  of  the  incident  power.  This  fraction  is  given  by 


A 

p  =  _2_ 

iN  ' 2tt'  4n 


'4-  IkV 


(5-14) 


-1  a  i 

where  8  =  2  tan  (t^— )  and  K  is  the  strength  of  the  magnetic  line 

s 

source.  A  normally  incident  plane  wave  in  region  a  produces  an 
intercepted  time  average  power  given  by  (|E*|  =1  volt/meter) 


P,M  =  w  n  |  H1 1  2 

iN  a  a  z 


(5-15) 


To  compute  the  time  average  power  into  region  c,  the  real  part  of  the 
Poynting  vector  flux  at  is  integrated  along  This  gives  the 

formula 


2  +3  c  *  -*-3 

P  =  -  vfY  1  VJ 

t3  k  n  1  1 

o  o 


(5-16) 


This  is  the  usual  equation  for  power  flow  into  a  network  represented  by 

Q 

an  admittance  matrix  [H  ]  except  for  the  factor  of  1/k  n  •  This  factor 

o  o 

appears  because  of  the  factor  -  k  r|  in  Eq.  (2-19).  Another  method  for 

o  o 

finding  P^  is  to  integrate  the  Poynting  vector  flux  over  the  far  field 
in  region  c.  This  gives 


46 


Pt3  "  Re 


tt/2 

— tt/2 
tt/2 


EC  x  hC*  •  r  rd(J) 


n  2tt  J 


|ief  V3|2  d<{) 


-tt/2 


(5-17) 


w 

For  the  TM  case,  an  electric  line  source  is  placed  at  (xs>  -y) 
in  region  a.  The  fraction  of  the  time  average  power  intercepted  by  the 
slit  is  given  by 


kn  0k  n 

p  =  e-L)  ..,a_a  I  ,i  I  2  _ a_a 

riN  '2ti  4  1  1  8  it 


ri  i  2 


(5-18) 


where  I  is  the  strength  of  the  electric  line  source  and  0  is  the  same 
as  before.  A  normally  incident  plane  wave  gives 


_  |  i,2 

piN  n  |Ezl 

a 


n 


(5-19) 


The  time  average  power  P  for  the  TM  case  is  given  by  Eqs.  (5-16)  and 
(5-17)  but  with  TM  quantities  used. 


5.3.  Power  Gain  and  Normalized  Field  Pattern 


The  power  gain  function  is  defined  as  the  ratio  of  the  radiation 
intensity  in  a  given  direction  to  the  radiation  intensity  which  would 
exist  if  the  transmitted  power  were  to  radiate  uniformly  over  half  space. 
This  gives,  for  the  TE  case 


G(<())  =  lim  - 
r  ->  oo 
m 


TTr  n  |  Hc  I  2 
me1  z  1 


2n  P 

C  t 


|ief  v3|2 


and  for  the  TM  case, 


(5-20) 
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(5-21) 


G(<J>)  -  lim 


Trr  o 

_ ! 5.  |EC|2 
P_  1  z1 
m  c  t 


2n  P 
c  t 


I imf  v3!2 


If  the  medium  in  region  c  is  lossy,  the  term  "radiation  intensity" 
does  not  apply.  Here  it  is  desirable  to  define  a  directive  gain  as 
the  ratio  of  the  time  average  power  density  produced  at  a  certain 
range  in  a  certain  direction  to  the  uniformly  radiated  time  average 
power  density  at  that  range. 

The  normalized  field  pattern  is  another  quantity  of  interest 
and  is  computed  in  the  far  field  by 


Fe(<S>) 


for  the  TE  case,  and  by 


Fm  (4>) 


lHz (*) 1 

|fef  v3| 

IhcI 

z  max 

|ief  v3| 

'max 

|E>)|  _ 

|^nf  -3| 

|ec| 

1  z 1  max 

|?*f  v3| 

'max 

(5-22) 


(5-23) 


for  the  TM  case.  The  denominator  denotes  the  maximum  value  of  field 
magnitude  measured  in  region  c. 


5.4.  Slit  Impedance  and  Polarizability 

3 

Since  the  fields  in  region  c  arise  from  a  magnetic  current  M 
radiating  in  the  presence  of  a  complete  conductor  at  x  =  d,  the  magnetic 
vector  potential  for  these  fields  may  be  written  as 
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(5-24) 


F 


M3(y’)  H<2)(kjr  -  y’£|)dy’ 


where  the  coordinates  of  Figures  7  and  8  are  used  and  r  replaces  r  as 

—  m 

the  field  point.  The  following  assumptions  are  made: 


a)  w  =  w  =  w. 

a  c 

b)  kcy'  «  1  for  y'  £  [-  j]  . 

c)  The  field  point  r_  is  such  that  |rj  »  y'. 

d)  The  incident  field  is  constant  on 


The  Hankel  function  in  Eq.  (5-24)  is  now  approximated  by  [24,  Appendix  D] 

^kV  (1  +  1  Sin  *)(1  +  Jky'sin  ^)e  C  (5-25) 

2 

Neglecting  terms  of  order  y'  and  higher,  one  can  write  the  magnetic 
vector  potential  as 
-jk  r  w 


F  - 


/2irjk  r 

c  w 


2  M3(y,)dy’ 


-jk  r 

jk  e  sin  <f> 


/2rrJiTr 


w 

2 


2  M3(y*)y’(l  +  2jT7)dy’ 


(5-26) 


The  first  term  in  the  above  equation  is  dominant,  and  hence  will  be  con¬ 
sidered  in  the  following  discussion  for  each  polarization. 
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For  the  TE  case,  the  magnetic  field  in  region  c  is  given  by 


jk 


HC  = - 

z  n  z 

c 


-  F  =  -  — 


jk  -jk  r 
_ c  e  c 

n  \J  2irr 


w 

2  M2(y')dy' 


J 

w 

2 


A  slit  impedance  Z  is  defined  by  the  equation 


w 

j2  M2(y’)dy'=  Z(H^C  w) 
w 

“  2 

isc 

where  H  ^  is  the  short  circuit  incident  magnetic  field  evaluated  at 
the  center  of  .  This  is  obtained  by  letting  the  sources  in  region 
a  radiate  in  the  presence  of  a  complete  conducting  plane  at  x  =  0. 

Q 

The  magnetic  field  is  now  written  in  the  form 


HC  = 

2  nc 


jk  -jk  r 
c  c  „isc 

2n  e  HZ1 


which  may  be  compared  to  the  measurement  vector  formula  given  in 
Eq.  (5-3). 


For  the  TM  case,  the  electric  field  in  region  c  is  given  by 


EC  = 
z 


-  z 


V  x  F  ~  -  cos  4> 


■^'jN 


(y’)dy' 


If  is  the  magnetic  charge  density  associated  with  the  transverse 
directed  magnetic  current  then 


(5-27) 


(5-28) 


(5-29) 


(5-30) 
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2  My(y')dy'  = 
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M3)dy ' 
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w 
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=  jto 


2  ,  ,  . 
y  P  dy' 
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w 

2 


=  jwy  p 

c  my 


(5-31) 


where  the  continuity  equation 


V'  •  MJ  =  -  jtop 
—  —  m 

has  been  used  and  a  magnetic  dipole  moment  is  defined  by 


*m 


jwu 


w 

ry  o 

r(y')dy' 


w 

2 


(5-32) 


(5-33) 


A  slit  polarizability  is  now  defined  in  terms  of  the  y  component  of 
Eq.  (5-33)  by  the  equation 


p  =  a  (wH1®0) 
my  m  yl 


(5-34) 


where  again  is  the  tangential  component  of  the  short  circuit  inci¬ 

dent  magnetic  field  evaluated  at  the  center  of  T  .  The  electric  field 
is  expressed  in  terms  of  this  polarizability  as 


E 


c 

z 


k  »n 

c  c 


2 


cos  <J> 


(5-35) 


which  also  may  be  compared  to  the  measurement  vector  formula  given  by 

Eq.  (5-9).  In  the  preceeding  discussion,  Z  has  the  dimension  of  ohms 

and  oi  has  the  dimension  of  length, 
m 
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Chapter  6 


MODAL  SOLUTION 

6.1.  Problem  Formulation 

The  discussion  thus  far  has  been  concerned  with  slits  of  arbi¬ 
trary  cross  section.  If  the  slit  has  a  rectangular  cross  section,  the 
equivalence  principle  may  be  used  to  break  the  problem  up  into  three 
isolated  regions  [16].  The  fields  in  region  b  then  have  a  modal  repre¬ 
sentation,  which  is  not  necessarily  the  case  when  the  slit  cross  section 
is  an  arbitrary  shape. 

In  this  chapter,  a  special  case  of  slit  cross  section  is  con¬ 
sidered  which  allows  region  b  to  be  considered  as  a  sequence  of  two- 
dimensional  rectangular  cavities  each  coupled  to  another  by  an  aperture. 
One  such  possible  configuration  is  shown  in  Fig.  9.  Each  cavity  may 
be  filled  with  different  material.  In  general,  region  b  will  consist  of 
Q  such  rectangular  sub-regions.  The  q  rectangular  sub-region  is 
specified  by  dimensions  as  shown  in  Fig.  10.  The  equivalence  principle 
is  applied  as  in  [16],  so  that  each  aperture  face  T  ,  for  q  =  1,2,..., 

Q  +  1,  is  replaced  by  a  perfect  electric  conductor  with  a  magnetic  current 
sheet  residing  on  each  side.  These  magnetic  current  sheets  are  chosen 
to  assure  continuity  of  the  tangential  component  of  electric  field  at 
each  r  in  the  original  problem. 

The  original  problem  has  now  been  broken  up  into  Q  +  2  isolated 
regions,  each  with  postulated  equivalent  sources.  The  total  field  in 
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Fig.  10.  The  qth  rectangular  cavity' of  region  b. 


region  a  is  (E  +  JE  ,  H.  ■+-  _H  ) ,  where  (_E  ,  H  )  is  due  to  sources  (.J  ,  M 
radiating  in  front  of  a  complete  conductor  at  x  =  0  and  (E3,  H3)  is  due 
to  magnetic  current  sheet  M  at  F^  radiating  in  front  of  a  complete  con¬ 
ductor  at  x  =  0.  The  magnetic  current  sheet  is  given  by 

Mx  =  x  x  E3  at  T  ( 

The  total  field  in  each  sub-region  bq  of  region  b  is  (E^ ,  H^)  due  to 

magnetic  currents  -  M  at  T  and  M  ,  at  F  ,,  radiating  inside  a  closed, 

-q  q  -q+1  q+1  b 

two-dimensional,  conducting  box  filled  with  material  of  permeability  u, 

bq 

and  permittivity  E^q*  These  magnetic  currents  are  given  by 


M  =  x  x  E 
-q  -  - 


M  ,  =  x  x  E  at  T 

—q+1  —  —  q+1 

c  c 

The  total  field  in  region  c  is  (E  ,  .H  )  due  to  the  equivalent  magnetic 
current  sheet  at  radiating  in  front  of  a  complete  conductor 

at  x  =  d.  This  magnetic  current  sheet  is  given  by 


Vi = 


The  preceding  equivalent  magnetic  current  sheets  must  also  be 
determined  so  as  to  assure  continuity  of  the  tangential  component  of 


magnetic  field  at  each  aperture  face  in  the  original  problem.  Thus  the 


following  equations  are  written: 


•  -  “t1! 


'•“’V  +  lJ5'2(«2)  +  ^-2(«2)  -  ^•3(M3)  -  0 


(6-5) 


-^lh’Q-S-l)  +  ^Qlb’Q(V  +  ■  ^Q,<H1<Vl>  ■  at  r 


+ ^r<v> +  2^i®’  Vi>  ■ a  «  fQ+i 

hs  tic  i 

The  operators  L  and  L  are  defined  by  Eq.  (2-2)  and  H  is  the 
— t  — t  — t 

tangential  component  of  the  magnetic  field  due  to  sources  (J1,  M1) 

radiating  in  the  half  space  filled  with  (y  ,  £  )  in  the  presence  of  a 

a  a 

complete  conductor  at  x  =  0.  The  factors  of  2  in  Eqs.  (6-5)  come  from 

tnh  ti 

the  images  of  the  magnetic  current  sources.  The  operator  ’  gives 
the  tangential  component  of  magnetic  field  at  due  to  a  magnetic  cur¬ 
rent  sheet  at  F  radiating  inside  the  two-dimensional  cavity  bm.  The 
n 

superscript  h  denotes  magnetic  field  operator. 

Using  the  method  of  moments  procedure,  one  can  reduce  Eqs,(f>~5)  to  a 


system  of  matrix  equations.  Each  aperture  face  F^  is  divided  into 
segments  of  constant  length  A^.  It  is  assumed  that  the  equivalent  magnetic 
current  sheets  may  be  expanded  as  a  linear  combination  of  basis  functions 
defined  on  each  domain  T  .  This  is  written  as 

q 
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(6-6) 


M  = 

~q 


N 

f 

n=l 


V 

q,n 


e 

-q,n 


(y) 


on  r 

q 


for  q  =  1,2,...,  (W-l.  The  coefficients  V  are  complex  scalars  to 

q,n 

be  determined  and  the  functions  e^  ^  will  be  specified  for  each  polari¬ 
sation.  The  testing  functions  are  chosen  to  be  the  same  as  those  used 
for  expansion  and  a  suitable  symmetric  product  is  defined  by 


<A  ,  B  > 

-q  -q 


A  •  B  dy 

J  -q  -q 
Q 

u  r 
q=l  q 


(6-7) 


where  A  and  B  are  vector  functions  defined  on  F  . 

-q  -q  q 

Equation  (6-6)  is  substituted  into  Eqs.  (6-5)  for  the  unknown 

magnetic  currents.  The  symmetric  product  of  the  qtb  equation  of  (6-5)  is 

then  taken  with  e  for  m=  1,2,...,  N  and  q  =  1,2,...,  Q+l.  The 
— q ,  m  q 

result  is  a  system  of  linear  equations  which  is  written  in  matrix  form  as 
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where  V  are  column  vectors  containing  the  coefficients  V  ,V  . V 

q  q ,  1  q ,  2  q,N 


(6-8) 
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The  matrix  in  Eq.  (6-8)  is  block  tridiagonal.  Each  submatrix  has 
elements  given  by  the  formulas: 


Ya  =  - 
mn 

M  <e.  .  2Lba(0,  e.  )> 

o  o  1  ,m  t  —  — l,n 

(6-9) 

(Ybq) 
''ll  mn 

=  -  k  n  <e  ,  Bqh,q (e  )> 

o  o  - q,m  — tq  — q ,n 

(6-10) 

(Ybq) 

V  12 'ran 

=  k  n  <e  ,  Bjh»q+1(e  )> 

o  o  ~q  ,m  — tq  — q+l,n 

(6-11) 

(Ybq) 

=  Vo  V,«’  ^(^,n}> 

(6-12) 

(Ybq) 

U22'mn 

=  -  Vo  V,m-  ^q;?+1(V 

(6-13) 

YC  =  - 
mn 

kono  <— Q+l,m’  2^tQ+l(-’  ^Q+l,nJ> 

(6-14) 

where  the  factor  -k  n  has  been  multiplied  through  Eqs.  (6-5).  The 
o  o 

non-zero  elements  of  the  excitation  column  are  given  by 


I1  =  k  n  <e  >  2h\  > 

m  o  o  — ]  ,m  — tl 


(6-15) 


The  ranges  of  the  indices  m  and  n  are  determined  after  the  specification 
of  the  expansion  functions. 

The  original  problem  may  be  viewed  as  Q+2  cascaded  networks.  Each 
network  represents  an  isolated  region  and  has  a  characteristic  admittance 
matrix  which  depends  only  upon  the  properties  of  the  region  which  the 
network  represents.  Due  to  the  choice  of  expansion  functions  made  in  the 
next  two  sections,  it  will  be  seen  that  the  sub-matrices  of  the  composite 
matrix  in  Eq.  (6-8)  have  the  following  properties: 
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1) 


3  C 

[Y  ]  and  [Y  ]  are  symmetric  Toeplitz  matrices,  hence 
only  one  column  of  each  need  be  computed. 

2)  Y^  and  Y^  ate  symmetric  matrices. 

3)  (vjj)  =  (Y?J)  . 

iz  mn  L i  nm 

The  matrix  of  coefficients  in  Eq.  (6-8)  is  block-tridiagonal  and 
symmetric,  hence  special  computational  and  storage  procedures  are 
used . 


6.2.  TE  Case 

For  this  polarization,  the  equivalent  magnetic  current  sheets 
are  z  directed.  Thus  expansion  functions  for  the  aperture  faces  T 

q 

and  T  , .  are  chosen  as 
q+1 


e  (y)  = 
— q,n 


for  n  =  1,2, ... ,  N  and 

q 


lz  (n-1)  A  <  y  -  y„  <  nA 

q  -  £q  -  q 


y  elsewhere 


flz  (n-1)  A  <  y  -  y  <  nA 

q+1-  rq  —  q+1 


Vi,»<y) 


(6- 


(6- 


y  elsewhere 


for  n  =  1,2, ... ,  N 


Also,  in  the  above,  A  =  w  /N  for  q  =  1,2,..., Q+1 


q+1'  '  '  q  q  q 

Equations  (6-16)  are  written  with  respect  to  the  coordinate  system  shown 

in  Fig.  10.  Figure  11  shows  the  order  in  which  the  functions  e  appear 

-q,n 

on  aperture  face 


16.1) 


16.2) 
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everywhere  in  region  bq  except  where  the  sources  are.  Considering  the 

magnetic  current  sources  e  at  T  and  e  ,,  at  T  ,  separately,  one 

— q,n  q  — q+l,n  q+1  r  1 ' 

may  solve  Eq.  (6-20)  with  the  appropriate  boundary  conditions  on  the 
cavity  walls.  The  resulting  magnetic  field  operators  are  then  de¬ 
fined  by 


Bqh’q(e  ) 
-z  ~ q,n 


jkb 


00  £  cos  k4  (x  -  dn  ) 

i  p  xP  q 


kq^bq  p=0  kq  sin  kq  d 
xp  xp  q 
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p^y 


and 
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e  _(y')  cos  dy' 

J  q.n  h 

n  q 


ml 


(6-21) 


Bqh',»+1(.  ,  ) 

~~ z  — q+1 ,  n 


=  z 


^kb 
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q  h_ 


£L_  Y  -2 _ ^ 


hqnbq  p=0  kq  sin  lei  dn  “q 
xp  *P  q  M 


HZl  Hv' 


Vi,n(y,)  cos-^dy' 
q 


(6-22) 


where  the  subscript  t  has  been  replaced  by  z.  The  second  subscript  on  B 

will  be  either  q(x  =  0)  or  q+l(x  =  d  ).  In  the  above,  £  =  1  for 

q  P 

p  =  0,  =  2  for  p  ^  0,  and 


(kq  )2  =  k? 

XP  DC 


<f>2 

q 


Thus,  for  Eqs.  (6-10)  and  (6-13),  one  may  obtain 


(6-23) 
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"  *"  '  • 


y  -  y. 


rq 


q+i 


+  l-n)£  (n-l)Aq+1  <  y  -  yrn  <  nA 


rq  -  q+1 


Vl.n(y) 


y  -  y 

/  rq 

Aq+i 


+  l+n)£ 


nAq+1  £  y  -  yrq  <  (n+l)Aq+1  (6-29.2) 


3 

r- 

l-’ 


0  y  elsewhere 

for  n  =  1,2,. . . ,  N  -1.  In  the  above,  A  =  w  /N  for  q  =  1,2,...,  N  -1. 

q+1  q  q  q  q 

Equations  (6-29)  are  written  with  respect  to  the  coordinate  system  shown 

in  Fig.  10.  Figure  12  shows  the  order  in  which  the  functions  e  appear 

-q,n 

on  aperture  face  T  . 

The  elements  of  [Y  ]  are  determined  from  a  specialization  of 
Eq.  (4-15)  to  ri.  This  is  written  as 


M. 


ra  o  o 


mn  k  ri 
a  a 


{S(m-1,  n-1,  1,  1)  +  S (m-1 ,  n,  1,  -1) 


+  S(m,  n-1,  -1,  1)  +  S (m,  n,  -1,  -1)} 


(6-30) 


kJL 


for  1  £  m  £  N^,  1  £  n  £  and  S  is  defined  by 

1  1 


S(m,n,p,q) 


(6-31) 


j  [  +  ¥  (^T  +  b  ~  ^Ho2)  <yI I  ~  T  +  (m-n)  I  >du’du 

if  m  ^  n 
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r* 


where  y  =  k  A,  and  the  function  a  is  defined  by  Eq.  (3-17).  The 

matrix  [Y  ]  is  determined  from  Eq.  (6-30)  where  y  =  k  A.,.,  k 

c  Q+l  c 

replaces  k  ,  and  n  replaces  ri  . 

a  c  a 

To  compute  Eqs.  (6-10)  -  (6-13),  TM  fields  must  be  constructed 
inside  a  two-dimensional  perfectly  conducting  cavity.  To  do  this,  an 
electric  vector  potential  Abq  =  \J;bqz  is  chosen  for  each  region  bq. 

The  fields  in  each  region  are  obtained  by 

l1”  '  -  A”’  «-32> 

Hbq  .  V  *  Abt>  -  S  ^  -  2  <6-33> 

Each  scalar  function  ipb<^  satisfies  Eq.  (6-20)  everywhere  in  region  bq 

except  where  the  sources  are.  Considering  the  magnetic  current  sources 

e  at  T  and  e  ,  at  T  ,  separately,  one  may  solve  Eq .  (6-20)  with 
— q , n  q  —q+l , n  q+l 

the  appropriate  boundary  conditions  on  the  cavity  walls.  The  resulting 
y  components  of  the  magnetic  field  operators  are  defined  by 


Bqh,q(e  )  =  l 
-y  -q,n  JL 


JLL 
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(6-34) 
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eq+l,n(y,)  sin-^dy- 


(6-35) 
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where  the  subscript  t  has  been  replaced  by  y.  The  second  subscript  will 

be  either  q(x=0)  or  q+l(x=d  ).  Equation  (6-23)  also  holds  for  k  .  Thu 

q  xp 

for  Eqs.  (6-10)  and  (6-13),  one  may  obtain 


u  2  jo  A  °°  / - — 

(Y.  ?)  =  — — — -  (k  A  )— 3-  £  1  -  (  -^  )  cot  kq  d 

11  mn  nbq  o  q  hq  kbqhq  xp  q 


Sln  T  (y£q  +  “V 
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*  Sln  T  (y£q  +  nAq) 
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for  1  <  m  <  N  -1,  1  <  n  <  N  -1,  and 
-  -  q  ~  ~  q 


2  i  ri  A  00  i — 
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q  q 


for  1  .<  m  Nq+^-l,  1  n  <_  Nq+^-l.  The  formula  for  Eq.  (6-11)  is 
given  by 


(YU>„n  '  nbq°  'VW  X/1  '  (i§r>T  CSC  k’pdq 


a+1  •,  2 


Sin  h“  (y*q  *  "V 


sin  f—  (y  +  nA  ) 
hq  rq  q+1 


(6-38) 


for  1  £  m  _<  Nq-1,  1  £  n  N  ^-1.  The  formula  for  Eq.  (6-12)  is  easily 
seen  to  be 


(Y«?)  =  (Y*J) 

21  mn-  12  nm 


(6-39) 


for  1  <  m  <  N  ,.-1,  1  <  n  <  N  -1. 

-  -  q+i  -  -  q 

The  non-zero  elements  of  the  excitation  column  vector  are  ob¬ 
tained  as  in  Section  4.3.  Thus  for  an  incident  plane  wave,  Eq.  (4-33) 


becomes 


1  i 

.  2k  n  jk,K  -  sin  <J) 
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(6-40) 


A  formula  for  line  source  incidence  may  be  obtained  by  a  specializa¬ 


tion  of  Eq.  (4-36). 


It 


Chapter  7 


APPROXIMATE  SOLUTION  TO  SLIT  OF  RECTANGULAR 
CROSS  SECTION 


7.1.  Problem  Specialization 


In  this  chapter,  a  slit  of  rectangular  cross  section  is  con¬ 
sidered  as  shown  in  Fig.  13.  This  is  a  special  case  of  the  more  general 
problem  considered  in  Chapter  6  and  Eqs.  (6-8)  become 
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(7-1) 


since  region  b  is  represented  by  one  rectangular  region.  This  region 
has  the  composite  admittance  matrix  representation  [Yb]  denoted  by 
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Equations  (7-3)  and  (7-4)  come  from  the  specialization  of  Eqs.  (6-10)  and 
(6-12)  and  the  last  two  come  from  the  fact  that  each  aperture  face  is 
divided  into  N  subsections  each  of  length  A  =  w/n.  The  elements  of  [Y3] , 
[YC],  and  I1  in  Eq.  (7-1)  are  obtained  from  Eqs.  (6-9),  (6-14),  and  (6-15) 
respectively. 

It  is  desired,  in  this  chapter,  to  obtain  a  suitable  approximation 
to  the  magnetic  field  inside  region  b  so  that  simplified  formulas  for 
Eqs.  (7-3)  -  (7-6)  may  be  obtained.  For  each  polarization,  certain  assump¬ 
tions  are  considered  so  that  the  fields  in  region  b  are  approximated  by 
plane  waves  travelling  in  the  +x  and  -x  direction.  This  corresponds  to 
neglecting  a  transverse  component  of  field  in  each  case  and  the  approxi¬ 
mation  is  expected  to  work  best  when: 

1)  The  incident  field  is  a  normally  incident  plane  wave. 

2)  Region  b  is  filled  with  a  dense  material,  i.e.,  the 
intrinsic  wavelength  of  the  material  filling  the  slit 
is  much  less  than  that  of  free  space. 

1.2.  TE  Case 

For  the  TE  case  the  magnetic  current  sheets  and  are  z  directed. 
If  they  are  constant  on  F^  and  respectively,  then  from  image  theory, 
the  fields  in  region  b  are  given  by  and  -x  travelling  plane  waves.  The 
currents  and  become  unbounded  at  the  ends  of  and  but  this  is  a 
very  localized  phenomenon,  particularly  when  region  b  is  a  dense  material. 
Thus  if  w  >>  X^,  the  magnetic  currents  should  be  constant  enough  to  neglect 
the  x  component  of  electric  field. 
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7.3.  TM  Case 


For  this  polarization,  the  magnetic  current  sheets  are  y  directed 
and  even  if  and  are  assumed  constant  over  and  the  field  in¬ 
side  region  b  is  not  represented  exactly  by  a  plane  wave.  If  one  con¬ 
siders  the  constant  magnetic  current  distribution  =  K  £  radiating 
inside  the  closed  conducting  box  0_<x<^d,  0  <  y  <  v,  then  one  may 
express  the  fields  inside  the  box  as 


Sln  V1  -  &  V*-d)) 


Sln  y1  ' 
sin  (  fl 


)  k,  (x— d  )  ) 


sin.  / 1 


V1  -  V 


Since  has  no  average  value  on  0  y  w,  and  if  k^w  >>  1,  it  is 

reasonable  that  one  could  neglect  as  compared  to  E^  and  H^. 

If  only  E  and  H  are  considered  in  region  b,  then,  one  may  pro 
z  y 

ceed  with  the  plane  wave  assumption  as  before  to  obtain 


1h  1  e,  (y)  C0S 

’  (e,)  =  j  — - -  (tv-1)  <  y  <  (n+l)A 

—y  — In  fl^  sm  k^d  —  — 

where  given  by  Eqs.  (6-29).  Thus  for  Eqs.  (7-3)  and  (7-4), 

one  obtains 
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where 


if  m  -  n  =  1 


if  m  =  n 


otherwise 


(7-17) 


Alternatively,  since  E  =0,  and  E  =  0  at  y  =  (n-l)A  and 

x  z 

y  =  (n+l)A  for  each  expansion  function  .e  ,  conducting  planes  may  be 

placed  at  the  ends  of  each  e,  so  that  each  e,  radiates  inside  the 

—In  —In 

closed  conducting  box  (n-l)A _<  y  <  (n+l)A.  Specializing  Eqs.  (6-34) 
and  (6-35)  to  this  problem,  one  obtains 


tu  i  -o  00  k  cos  k  (x-d) 
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(7-18) 


where  k^  =  k^  -  (^y)^-  Keeping  only  the  p  =  1  term,  one  obtains  for 
Eqs.  (7-3)  and  (7-4)  the  following  formulas 
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Chapter  8 


NUMERICAL  RESULTS  AND  DISCUSSION 

8.1.  Modal  and  Non-Modal  Results 

Some  representative  computations  are  presented  here  for 
filled  slits  of  various  cross  sections.  The  cross  section  shapes 
and  their  dimensions  are  shown  in  Figs.  14-16.  A  comparison  will 
first  be  made  between  the  results  obtained  from  the  non-modal  solu¬ 
tion  of  Chapters  1  to  4  and  the  modal  solution  of  Chapter  6.  The 
quantities  to  be  compared  are  the  magnetic  current  on  the  aper¬ 
ture  faces  ri  and  F^  and  the  electric  current  on  the  contour  C 
defining  the  slit  cross  section.  This  electric  current,  which  appears 
as  an  unknown  in  the  non-modal  solution,  is  computed  in  the  modal 
solution  by  first  solving  for  the  magnetic  currents.  Then,  using 
parallel  plate  guide  modes,  the  tangential  magnetic  field  due  to 
these  magnetic  currents  radiating  inside  a  closed  two-dimensional 
conducting  box  is  computed  just  inside  the  contour  C.  This  corre¬ 
sponds  directly  to  the  electric  current  unknown  in  the  non-modal 

solution.  The  magnitudes  of  the  electric  current  times  T)  and  the 

o 

magnetic  current  are  plotted  as  a  function  of  their  position  on  the 
contour  C  (cf.  Fig.  4).  Horizontal  and  vertical  lines  are  used  to 
denote  quantities  expanded  by  pulses  while  straight  lines  connect  the 
points  given  by  the  triangle  expansion  functions.  Each  line  F^  com¬ 
prising  the  contour  C  is  broken  up  into  N^+^  -  N^  straight  line  seg¬ 
ments  which  are  shown  scaled  in  the  figures.  The  excitation  consists 


77 


of  a  plane  wave  incident  from  the  left  of  the  slits  at  an  angle  4>i 

measured  with  respect  to  the  negative  x  axis. 

The  first  example  to  be  considered  is  a  slit  of  rectangular 

cross  section  of  width  w  =  0.8X  and  thickness  d  =  0.2X  (X  being  the 

o  o  o 

wavelength  in  free  space).  Regions  a,  b,  and  c  are  always  assumed  to 
be  free  space  unless  otherwise  specified.  Excellent  agreement  is  noted 
between  the  two  solutions  for  the  TE  case  when  e,  =  C  in  Figs.  17  and 

D  O 

18  and  when  e,  =  5e  in  Figs.  19  and  20.  A  numerical  check  is  made  on 
d  o 

the  solution,  if  region  b  is  lossless,  by  testing  to  see  if  conserva¬ 
tion  of  power  flow  holds  at  the  two  aperture  faces.  Thus  the  real 
parts  of  the  Poynting  vector  flux  at  F^,  denoted  by  P^,  and  ,  de¬ 
noted  by  P  ,  are  computed.  They  are  compared  by  the  number  AP  defined 
as 


AP  =  100 


P  -  P 
tl  t3 


(Ptl  +  Pt3>/2 


(8-1) 


which  should  be  as  small  as  possible.  This  number  is  computed  readily 
in  the  non-modal  solution  and  is  given,  when  appropriate,  in  the  figure 
captions  of  the  computed  examples.  Cases  for  which  the  computation 
(8-1)  is  not  reliable  is  when  the  slit  is  cut  off  in  the  TM  case  and 
the  imaginary  part  of  the  Poynting  vector  flux  is  much  greater  than  the 
real  part. 

Two  TM  cases  of  the  previous  example  are  shown  in  Figs.  21  and  22 
where  agreement  between  the  magnetic  currents  is  again  excellent.  The 
noticeable  difference  between  the  two  solutions  occurs  in  the  electric 
current  at  the  ends  of  the  illuminated  aperture  face  An  analysis 


of  the  field  near  the  edges  following  the  procedure  of  Meixner  [28] 
shows  that  the  magnetic  current  should  approach  zero  with  an  infinite 
slope  right  at  the  ends  of  Also,  the  electric  current  should 

become  singular  at  these  points.  An  explanation  as  to  why  the  com¬ 
puted  non-modal  TM  electric  current  doesn't  behave  as  expected  at 
the  ends  of  the  aperture  face  lies  in  the  fact  that  the  edge  behavior 
of  the  magnetic  current  is  not  represented  adequately  by  the  tri¬ 
angle  expansion  functions.  They  constrain  the  magnetic  current  to 
approach  zero  linearly  at  the  ends  of  the  aperture.  The  modal  compu¬ 
tation  of  this  electric  current,  however,  always  gives  a  zero  value 
of  electric  current  at  y  =  0  and  y  =  w  because  of  the  nature  of  the 
modal  expansion.  This  modal  computation  converges  non-uniformly  to 
the  electric  current  at  these  points. 

To  illustrate  a  more  detailed  behavior  of  the  aperture  fields, 

a  smaller  slit  (w  =  0.2X  ,  d  =  0.2A  )  is  considered  where  the  number 

o  o 

of  unknowns  is  the  same  as  in  the  previous  example.  In  Figs.  23-26, 
for  the  TE  case,  the  agreement  in  the  electric  current  from  the  two 
solutions  is  still  excellent.  For  the  magnetic  current,  however,  a 
slight  oscillation  is  present  in  both  solutions,  the  greater  vari¬ 
ation  being  noticed  in  the  non-modal  results.  For  a  slit  of  this 
size,  it  is  not  expected  from  physical  considerations  that  the  true 
solution  would  have  any  oscillatory  behavior.  It  is  thus  concluded 
that  this  anomaly  is  a  characteristic  of  the  numerical  solution  even 
though  the  numerical  values  for  AP  are  small.  The  TM  cases  for  this 
example  are  shown  in  Figs.  27-29  where  the  modal  electric  current  was 
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csmputed  at  50  points  along  each  line  T  of  C.  Results  from  both 
of  the  solutions  compare  well  except  for  J*3  at  each  aperture  face 
and  r3.  The  oscillatory  behavior  in  the  electric  current  from  the 
modal  solution  is  due  to  the  triangular  representation  of  the  mag¬ 
netic  current  which  has  a  pulse  representation  of  magnetic  charge. 
This  results  in  a  logarithmic  singularity  in  the  y  component  of 
magnetic  field  at  each  break  point  of  the  triangles  and  hence 
accounts  for  the  difficulty  in  computing  the  electric  currents  on 
ri  and  r3  by  parallel  plate  guide  modes.  The  agreement  between 
electric  currents  on  V ^  and  is  quite  good  and  the  expected 
behavior  at  the  edges  is  observed.  For  the  x  component  of  magnetic 
field,  the  modal  series  converges  quite  rapidly.  The  effects  of 
increasing  the  dielectric  constant  in  region  b  until  there  is  a 
propagating  TM  mode  are  also  noticed  in  Figs.  27-29. 

The  preceding  results  for  slits  of  rectangular  cross  sections, 
which  have  received  extensive  consideration  in  the  literature,  gives 
one  an  idea  of  how  these  two  different  solutions  compare  for  this 
simple  cross  section.  As  a  further  check,  the  modal  and  non-modal 
solutions  for  the  slit  cross  section  of  Fig.  15  are  compared.  Here 
we  have  a  cross  section  composed  of  two  rectangular  regions  which 
may  be  solved  readily  by  either  method.  For  the  TE  case  shown  in 
Figs.  30  and  31,  excellent  agreement  is  obtained  for  all  quantities 
computed.  Slight  oscillations  are  again  noticed  in  the  results  for 
the  magnetic  currents.  Some  loss  is  added-  to  the  material  filling 
region  b  (e^  =  (1— j ) )  in  Figs.  32  and  33  which  corresponds  to  a 
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conductivity  of  0.017  siemens/meter.  Agreement  between  the  two 
solutions  is  again  quite  good.  The  TM  results  for  this  example 
are  shown  in  Figs.  34-37  where  the  magnetic  currents  from  the 
two  solutions  are  seen  to  agree  quite  closely.  There  is  more 
discrepancy  between  the  electric  currents  where  the 
largest  differences  occur  right  at  the  corner  of  the  cross 
section  on  line  T To  improve  this,  more  expansion  functions 
are  needed  on  since  the  introduction  of  an  edge  there  causes 
tremendous  changes  in  the  fields  close  to  the  edge. 

Allowing  for  a  few  discrepancies  at  some  points  on  the 
contour  C,  agreement  between  the  computed  quantities  in  the  pre¬ 
vious  examples  is,  in  general,  quite  good.  In  fact,  as  far  as 
the  transmitted  fields  in  region  c  are  concerned,  there  is  negli¬ 
gible  difference  between  the  two  solutions  for  a  given  problem. 
Figures  38-41  show  some  plots  of  transmission  coefficient  versus 
plane  wave  angle  of  incidence  (as  measured  from  the  negative  x-axis) 
for  some  of  the  previously  considered  slits. 

The  primary  advantage  of  the  non-modal  solution  is  that  it 
may  be  used  to  solve  slits  whose  cross  sections  cannot  be  broken  up 
exactly  into  rectangular  sub-regions.  Figures  16a-h  show  some  slits 
of  various  cross  sections,  some  of  which  are  similar.  Figures  42 
and  43  show  the  effects  on  the  radiation  patterns  in  region  c  if 
the  cross  section  is  changed  so  that  some  of  the  conducting  plane 
pushes  into  the  rectangular  cross  section.  These  perturbed  cross 
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sections  are  shown  in  Figs.  16b, c,  and  e.  The  overall  effect  is 
to  shift  the  patterns  slightly  upward  in  the  TE  case  while  almost 
no  effect  is  noticed  in  the  TM  case.  If  the  rectangular  cross 
section  is  changed  so  that  some  of  the  conducting  plane  material 
is  carved  out,  thus  making  region  b  bigger,  a  much  more  noticeable 
upward  shifting  of  the  pattern  is  noticed  in  Fig.  44.  There  is 
again  a  negligible  change  in  the  TM  slit  patterns  for  those  cases. 
Lastly,  Fig.  45  shows  the  patterns  for  the  slits  of  Figs.  16d ,  f 
and  g  when  the  plane  wave  is  incident  at  an  angle  of  45°. 
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Fig  •  17  •  TE  Plane  wave  incident  at  (J)  =  0°  on  rectangular 
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Fig.  18.  TE  Plane  wave  incident  at  4>  =  45  on  rectangular  slit 

d  =  ,2A  ,  e  =  e  ,  AP  =  0.01%  ,  J  =  n  Jb. 
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Fig.  19  .  TE  plane  wave  incident  at  cj>  =0°  on  rectangular  slit 
d  =  .2X  ,  eL  =  5e  ,  AP  =  0.50%  ,  J  =  n  Jb. 


Fig.  21.  TM  plane  wave  incident  at  <f>  =0°  on  rectangular  slit, 


3.0 
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Fig.  22.  TM  plane  wave  incident  at  $  =  45°  on  rectangular  slit, 

d  =  .2X  ,  e.  =  5e  ,  AP  =  0.02%,  J  =  n  Jb. 
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Fig.  24.  TE  plane  wave  incident  at  <p  =  45°  on  rectangular  slit, 
d  =  . 2A  ,  =  e  ,  AP  =  0.35%  ,  J  =  n  A 
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Fig.  25.  TE  plane  wave  i  ident  at  <J>  =  0°  on  rectangular  slit 

d  =  .2X  ,  e,  =  ,  AP  =  3.10%  ,  J  =  n  Jb- 


0*9 


Fig.  26.  TE  plane  wave  incident  at  <p  =  45°  on  rectangular 


Fig.  28.  TM  plane  wave  incident  at  4>  =  0°  on  rectangular 


wave  incident  at  4>  =0°  on  rectangular 


Fig.  31.  TE  plane  wave  incident  at 


Fig.  32.  TE  plane  wave  incident  at  <J>  .=  0°  on  slit  of  Fig.  15,  e  =  (l-j)e 
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Fig.  2  .  TM  plane  wave  incident  at  (f>  =  60°  on  slit  of  Fig.  15,  e  =  c  ,  AP  =  17.25%  , 
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Fig.  38.  TE  Transmission  coefficients  vs.  (j)1  for  slit  of 
Fig.  15  (circles)  and  rectangular  slit  w  =  .  8Aq, 

d  =  .5X  ,£,*£. 

o  b  o 


Fig.  39.  TE  transmission  coefficients  vs.  4)1  for  rectangular 
slit  w  =  .2\  ,  d  =  .2A  (circles)  and  rectangular 


1 
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Fig.  40.  TM  transmission  coefficient  vs.  ij)1  for  slit  of 
Fig.  15,  eb  =  eo. 
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Fig.  41.  TM  transmission  coefficient  vs.  <|>  for  rectangular 

slit  w  =  .2X  ,  d  =  .2X  ,£,=£. 

o  o  b  o 


[I 


18 


8.2.  Slit  Impedance  and  Equivalent  Circuit  for  TE  Case 


When  the  slit  cross  section  is  rectangular,  several  approximate 
solutions  for  the  transmission  characteristics  are  readily  available 
after  making  suitable  assumptions.  The  first  cases  considered  are 
when  regions  a,  b,  and  c  are  filled  with  free  space  and  the  slit  width 
w  is  such  that  kQw  «  1.  As  discussed  in  Chapter  5,  the  transmission 
characteristics  may  then  be  determined  from  a  slit  impedance  for  the 
TE  case,  and  a  slit  polarizability  for  the  TM  case.  The  TE  case  is 
very  important  in  electromagnetic  shielding  problems  since  there  is 
always  a  propagating  mode  in  the  waveguide  region  b.  If  the  condi¬ 
tion  k^w  «  1  is  satisfied  and  d  >  w  then  the  two  assumptions  are  made: 

1)  Each  half  space  may  be  represented  by  a  lumped 
impedance. 

2)  Region  b  behaves  like  a  uniform  transmission  line. 

This  suggests  the  equivalent  circuit  shown  in  Fig.  46.  Regions  a  and  c 


Fig.  46.  Equivalent  circuit  model  for  rectangular  slit  when  k  w  <<1. 
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are  identical  and  are  represented  by  the  lumped  impedance  Z  which 

a 

is  defined  by 


(8-2) 


The  excitation,  I  ,  is  equal  to  the  negative  of  the  electric  current 

SC 

which  would  exist  at  if  r.  were  completed  by  a  perfect  electric  con¬ 
ductor.  This  is  a  form  of  the  equivalence  principle  used  by 
Collin  [30].  The  voltages  and  are  obtained  by  integrating 
the  electric  ixeld  vmagneuic  currents;  across  the  aperture  faces 
Fj_  and  T  respectively.  It  is  desired  to  find  the  slit  impedance 
defined  by  Eq.  (5-28)  which  may  be  rewritten  in  terms  of  the  equi-  * 
valent  circuit  parameters  as 

V„ 


Z  = 


(8-3) 


sc 


From  transmission  line  theory,  one  may  easily  obtain 


Z_ 

n 


2z  '''cos  kd  +  j  (1+z  2) sin  fcd 


(8-4) 


where  z  -  Z  / n  and  n  is  the  characteristic  impedance  of  free  space, 
a  o  o 

The  circuit  parameter  Z  is  found  by  computing  the  impedance  for  an 

a 

infinite  flange  (d  °°)  when  k  w  «  1.  If  one  assumes  a  constant 

o 

electric  field  in  the  aperture  equal  to  unity  [24,  p.  180],  the  re¬ 
sult  is 


=  ykw 

o  w[tt  -  2 j  In  ( — y-)] 
2e 


(8-5) 


where  y  =  1.781  and  e  =  2.71828...  Another  formula,  based  on  a  static 
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solution  of  the  flange  problem,  is  given  by  [31] 


Z 

a 


X 

o 


yk  w 

w[tt  -  2 j  In  (—£-)] 


(8-6) 


For  comparison  purposes,  one  last  formula  is  obtained  which  uses 
the  assumption  that  the  electric  field  in  the  aperture  is  that 
obtained  for  a  slit  in  a  zero  thickness  conducting  plane  [32]  and 
is  given  by 


Z 

a 


X 

o 


yk  w 

w[tt  -  2j  In  (•—-)] 


(8-7) 


For  kQW  small  enough,  there  is  little  difference  between  these 
formulas. 

The  transmission  coefficient  may  also  be  obtained  from  the 
equivalent  circuit  by  determining  the  power  delivered  to  the  load  Z 

a 

at  the  distance  d  and  dividing  by  the  power  delivered  to  the  cir¬ 
cuit  by  I  when  Z  =  r|  .  Thus  one  may  obtain 
sc  a  o 

T  =  ~  | Z | 2  Re{l/Z  }  (8-8) 

1 1  d 

O 

Some  plots  of  Z / n  and  T  versus  d/XQ  are  shown  in  Figs.  47-49 
for  rectangular  slits  of  various  widths.  The  solid  curves  repre¬ 
sent  computations  from  the  equivalent  circuit  in  Fig.  46  where  Z  is 
determined  from  Eq.  (8-7).  Circles  and  triangles  are  used  to  mark 
results  obtained  from  the  modal  solution  of  the  rectangular  slit 
computed  for  each  value  of  w  and  d.  Agreement  is  excellent  for  slits 
whose  widths  are  even  as  large  as  0.2Xo<  For  wider  slits,  the  slit 
impedance  Z  loses  its  usefulness  since  it  can  no  longer  accurately 
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determine  the  transmitted  fields. 

It  is  interesting  to  note  that,  as  w  becomes  smaller,  the 

transmission  coefficient  tends  to  a  peak  value  when  d  is  close  to 

a  multiple  of  A^/2.  These  peaks  approach  a  value  of  1/wtt  for  small 

w  which  can  be  shown  from  the  equivalent  circuit.  This  behavior  is 

predictable  in  that  a  resonance  condition  is  achieved  when  the 

reactance  of  Z  ,  at  the  distance  d,  cancels  that  of  Z  at  the  source 
a  a 

Igc*  Thus,  even  though  very  little  incident  energy  impinges  upon 
the  slit  for  small  w,  the  transmission  coefficient  may  be  quite  large 
for  some  thickness  d  when  compared  to  the  d  =  0  case. 
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Slit  impedance  and  transmission  coefficient  vs.  d/A 

o 

for  w  =  0.1Xq.  Circles  and  triangles  denote  modal 
results,  solid  lines  denote  equivalent  circuit  model 


For  wider  rectangular  slits,  the  plane  wave  approximation  of 
Chapter  7  is  useful.  Plots  of  transmission  coefficient  versus  plane 
wave  angle  of  incidence  are  computed  using  this  approximation  (denot¬ 
ed  by  solid  lines)  and  using  the  modal  solution  (denoted  by  circles) 
using  a  sufficient  number  of  unknowns  on  each  aperture  face.  A 
relatively  narrow  TE  slit  (w  =  1.4Xq)  is  shown  in  Figs.  50a  and  b 
where  increasing  the  thickness  from  0.2Xq  to  0.4Xq  improves  the  plane 
wave  solution  somewhat.  Even  closer  agreement  is  noted  as  the  width 
is  increased  to  3-OX^  as  seen  in  Figs.  50c  and  d.  When  the  slit  is 
wide,  the  plane  wave  approximation  is  better  justified  from  the 
analysis  given  in  Chapter  7.  This  is  clearly  seen  in  Figs.  51a  and  b 

where  £,  =  30e  so  that  the  slit  is  about  5  —■  times  its  free  space 
bo  2 

electrical  width.  A  TM  slit  of  width  w  =  1.4X  is  shown  in  Figs. 

o 

52a  and  b  where  it  is  again  observed  that  increasing  the  thickness 

d  improves  the  plane  wave  solution.  The  wider  slit  w  =  3X  is  shown 

o 

in  Figs.  52c  and  d  for  the  TM  case.  The  effects  of  adding  loss  to  the 
material  filling  region  b  are  shown  in  Figs.  53  and  54. 


0.4X  ;  c)  w  =  3X  ,  d  =  0.2X  ;  d)  w  =  3X  ,  d  =  0.4X  .  Circles  denote  modal  solution 
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(b) 

TE  transmission  coefficients  vs.  if)*  for  e,  =  30e  , 

b  o 

a)  w  =  1.4X  ,  d  =  0.2X  ;  b)  w  =  1.4X  ,  d  =  0.4X  . 

o  o  o  o 

Circles  denote  modal  solution  and  solid  lines  denot 
plane  wave  approximation. 
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53.  TE  transmission  coefficient  vs.  (p  for 
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54.  TE  transmission  coefficient  vs.  4'1  for 


Chapter  9 


CONCLUSIONS  AND  RECOMMENDATIONS 


A  non-modal  formulation  has  been  developed  for  studying  elec¬ 
tromagnetic  transmission  through  a  slit  of  arbitrary  cross  section 
cut  in  a  perfectly  conducting  ground  plane  of  finite  thickness.  The 
term  "arbitrary  cross  section"  means  that  the  two  slit  faces,  repre¬ 
sented  by  lines  and  are  coincident  with  the  planes  x  =  0  and 
x  =  d,  respectively,  and  that  the  rest  of  the  cross  section,  lines 
^  and  r4,  may  be  specified  in  an  arbitrary  manner.  After  some  modi¬ 
fication,  which  is  not  undertaken  in  this  report,  the  formulation  can 
also  be  applied  to  cases  where  the  material  filling  the  slit,  and 
hence  lines  and  F^,  may  protrude  out  into  the  half  space  regions 
a  and  c. 

A  modal  solution  is  presented  to  treat  slits  with  a  cross 
section  which  may  be  represented  by  a  chain  of  two-dimensional  rec¬ 
tangular  cavities  which  are  coupled  to  each  other  by  an  aperture. 

This  chain  coupling  approach  yields  a  block  tri-diagonal  matrix  which 
allows  special  computational  considerations.  In  dealing  with  a  cross 
section  composed  of  more  than  one  rectangular  sub-region,  the  same 
number  of  parallel  plate  guide  modes  is  taken  in  each  sub-region  to 
compute  the  fields.  This  number  is  more  or  less  arbitrarily  picked 
while  checking  to  see  that  the  real  power  flow  across  each  aperture 
face  is  essentially  the  same. 

The  above  two  formulations,  which  are  also  applicable  when  the 
material  filling  the  slit  is  lossy,  are  compared  to  one  another  in 
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various  examples.  The  advantage  of  the  non-modal  solution  lies  in 
the  arbitrary  nature  of  the  specification  of  the  slit  cross  section 
whereas  the  modal  solution  may  be  used  when  the  material  filling 
the  slit  has  certain  types  of  inhomogeneities.  The  computational 
effort  is  usually  greater  in  the  non-modal  solution,  especially  as 
the  slit  cross  secton  becomes  larger.  The  convergence  of  the  numeri¬ 
cal  solutions  is  only  investigated  in  the  sense  that  the  mean  differ¬ 
ence  between  the  two  solutions  is  observed  to  decrease  as  the  number 
of  unknowns  is  increased. 

A  brief  exposition  of  some  approximate  solutions  has  also  been 
presented.  It  is  well  known  that  a  small  hole  in  a  conducting  plane 
may  be  represented  in  terms  of  an  electric  and  magnetic  dipole  moment. 
For  the  two-dimensional  slit,  the  TE  case  is  of  more  importance  than 
the  TM  case,  because  there  is  then  always  a  propagating  mode  in  the 
parallel  plate  cross  section  no  matter  how  small  the  width.  Here,  a 
slit  impedance  is  defined  from  which  one  may  compute  the  trans¬ 
mission  characteristics  of  the  slit.  Next,  an  equivalent  circuit  is 
postulated,  based  on  a  transmission  line  model,  in  terms  of  the 
aperture  impedance  of  an  infinite  flange.  It  is  found  that  this 
equivalent  circuit  accurately  models  the  slit  for  widths  as  large  as 
0.2  wavelengths.  As  the  slit  width  becomes  very  small,  the  trans¬ 
mission  coefficient  becomes  very  large  for  conducting  planes  whose 
thickness  corresponds  to  the  distances  at  which  the  equivalent  circuit 
resonates.  This  circuit  model,  obtained  for  slits  of  a  rectangular 
cross  section,  could  be  extended  to  more  complicated  cross  sections 
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by  using  results  from  waveguide  junction  theory. 

For  slits  whose  cross  section  is  again  rectangular  but  whose 

widths  are  large  with  respect  to  the  wavelength  in  the  material 

filler,  a  so-called  plane  wave  approximation  can  be  used.  This 

greatly  simplifies  the  block  tri-diagonal  matrix  which  arises  in 

the  modal  solution  and  allows  very  rapid  computation.  It  is  found 

that  this  approximation  yields  accurate  results  when  k  w  »  1  and 

b 

is  slightly  better  in  the  TE  case.  This  approximate  solution  could 
also  be  applied  to  slits  which  are  filled  with  layers  of  dense 
material. 

Three  computer  programs  have  been  written  which  implement  the 
non-modal,  modal,  and  approximate  solutions  discussed  in  this  report 
They  are  listed  and  documented  in  [34], 
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Appendix 


NUMERICAL  APPROXIMATION  OF  INTEGRALS 


90  •  * 


In  Chapters  3  and  A,  several  integrals  appear  which  must  be 
computed  numerically.  These  integrals  are  defined  on  straight  line 
segments  in  the  x-y  plane.  The  Hankel  functions  and  are 

encountered  in  the  integrands  and,  since  they  are  computed  using 
polynomial  approximations  [29],  it  is  natural  to  choose  a  quadra¬ 
ture  formula  which  gives  the  highest  degree  of  precision  for  inte¬ 
grating  polynomials.  Thus  the  w$ll#kof-wiw  Gaussian  quatiraturt: 
formulae  [33]  are  used.  With  reference  to  Fig.  5,  consider  the 
integral 

m+1 

In  =  F(|r-r*|)dt'  (A-l) 

t 

n 

where  r  is  a  field  point  not  on  AC  ,  r'  is  a  source  point  on  AC  ,  and 
—  n  —  n 

F  is  a  well  defined  function  on  AC^.  Rewriting  r*  as 

r'  =  R  +  t'£ 

—  — n  — n 


where  R  is  a  vector  from  the  origin  to  the  midpoint  of  AC  ore  may 
— n  n 

rewrite  (A-l)  as 


AC 

1  --r- 

n  2 


AC 


F(|r  -  u’^  -  Rj  )du 


(A-2) 


AC 


where  t*  =  —rp-  u’.  This  is  now  in  the  form  to  approximate  by 


AC  Q  f  r\\  AC  tew 

j=i 


(A-3) 
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m  MMiim 


where  uj^  are  the  roots  of  the  Legendre  polynomial  of  degree  Q 


VU3W)) 


,(Q) 


and  Aj '  are  the  weighting  coefficients  determined  by 


,(Q)  . 


J  [1  -  (u<Q))2][P’(u<Q))]2 


Tables  of  uf^  and  f°r  various  values  of  Q  are  readily  available 

[33,  p.  337].  The  foregoing  rule  can  be  applied  in  succession  to 

multiple  integrals.  For  example,  if  I  is  well  defined  on  AC  ,  then 

n  m 

we  have  the  formula 


:m+l 

1  dt  'V 

n  'v 

t 

m 


AC  AC 
m  n 

4 


Q  Q 


l  l 

i=l  j=l 


F( 


fu^  t 


AC 
_ n 

2 


u«» 


t 

— n 


+  R 
— mn 


) 


(A-4) 
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